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INTRODUCTION 
Nature  of  the  Study 

The  control  of  a  neutron  chain  reaction  Is  probably  the 
most  important  aspect  of  atomic  energy.  The  ability  to  control 
a  nuclear  reactor  permits  the  usefulness  of  such  a  device.  As 
ia  well  known,  the  controllability  of  a  neutron  chain  reaction 
is  primarily  due  to  the  existence  of  delayed  neutrons,  i.e., 
neutrons  which  are  born  a  considerable  length  of  time  (relative 
to  neutron  lifetimes)  after  the  fission  process.  The  time  be- 
havior of  a  nuclear  reactor  must  be  well  understood  in  order 
that  the  very  concept  of  nuclear  energy  not  become  a  menace  to 
mankind.  The  study  of  reactor  kinetics  has  progressed  from 
early  in  the  Manhattan  Project  and  many  attempts  have  been  made 
to  understand  and  solve  reactor  kinetics  problems.  In  many  in- 
stances, in  order  to  solve  the  kinetics  problems,  each  reactor 
must  be  treated  individually  and  even  then  the  problems  may  be 
extremely  difficult  or  impossible. 

The  usefulness  of  any  theory  lies  in  its  ability  to  predict 
mathematically  the  characteristic  behavior  with  which  it  is  con- 
cerned. Therefore,  kinetics  theory  must  be  able  to  describe  the 
time  behavior  of  a  neutron  chain  reaction.  This  paper  investi- 
gates the  characteristics,  limitations,  and  utility  of  some  of 
the  presently  used  theories.  The  theory  utilized  in  this  paper 
is  very  similar  to  that  discussed  in  treatises  such  as  Glasstone 
and  Edlund  (8).  The  theory  was  applied  to  reactor  systems  for 
the  prediction  of  the  neutron  density  time  dependence  where  both 


step  and  finite  changes  of  reactivity  were  assumed.  This  paper 
will  demonstrate  the  agreement  found  between  theoretical  and  ex- 
perimental determination  of  the  reactivity  worth  of  reactor  con- 
trol rods,  measured  both  with  the  positive  period  method  and 
with  the  rod  drop  method. 

The  particular  reactor  of  interest  here  is  a  TRIGA  reactor 
(see  PLATE  I)  (designed  and  built  by  the  General  Atomic  Division 
of  General  Dynamics  Corp.).  Since  Kansas  State  University  will 
install  a  TRIGA  Mark  II  reactor  in  1961,  it  is  important  to  know 
beforehand  the  nature  of  the  kinetic  analysis  that  will  be  re- 
quired to  predict  the  behavior  of  such  a  reactor.  Because  of 
the  extreme  variation  in  complexity  of  the  equations  for  the 
various  kinetic  models,  it  is  important  to  realize  the  simplest 
model  which  will  suitably  describe  the  time  behavior  of  the  re- 
ac  tor  . 

Since  almost  all  solutions  of  the  various  kinetic  equations 
must  be  solved  with  the  aid  of  a  computer,  either  digital  or 
analog,  it  is  important  to  realize  an  approximate  solution  which 
will  accurately  predict  the  kinetic  behavior  of  a  reactor. 
Therefore,  solutions  of  the  kinetic  equations  using  both  three 
and  six  delayed  neutron  groups  were  developed  and  solved  with 
the  aid  of  the  Kansas  State  University  IBM-6£0  computer. 

Since  the  TRIGA  reactor  has  a  control  rod  drop  time  of 
approximately  0.35  seconds  from  "full  out"  to  "full  in",  the 
step  change  in  reactivity  assumption  is  not  really  Justifiable 
when  considering  neutron  flux  predictions  for  times  comparable 
to  the  rod  drop  time.  A  theoretical  development  is  given  for 


EXPLANATION  OF  PLATE  I 

General  view  of  the  TRIGA  reactor, 
looking  into  the  reactor  tank  at 
the  core.  The  ionization  chambers 
are  shown  just  outside  of  the  core 
reflector. 


PLATE  I 
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reactivity  insertion  as  a  function  of  time.  The  resulting  ana- 

lytical prediction  for  the  neutron  flux  as  a  function  of  time 

computed  from  the  one  group  model  was  compared  with  the  experi- 

mentally measured  neutron  flux  after  a  change  in  reactivity. 

Nomenclature 

r\ 

-  Neutron  density,  neutrons /cm'. 

£ 

-  Thermal  neutron  flux,  neutrons/cm  .  sec . 

D 

-  Diffusion  coefficient,  cm. 

V* 

-  Laplacian  operator 

r 

-  Spatial  coordinate 

t 

-  Time,  seconds 

f 

-  Fraction  of  total  neutrons  that  are  delayed 

e, 

-  Fraction  of  neutrons  in  the  ith  delayed  group 

Kh 

-  Effective  multiplication  constant 

te 

-  Infinite  medium  multiplication  constant 

*i 

-  Decay  constant  of  precursor  of  the  delayed  neutrons  of 

the  ith  group 

t, 

-  Precursor  density  of  ith  group,  atomic  nuclei/cm^. 

i 

-  Prompt  neutron  lifetime 

B* 

-  Buckling 

£* 

-  Macroscopic  absorption  cross  section,  cm"1. 

t 

-  Fermi  age,  cm  , 

V 

-  Neutron  velocity,  cm./sec. 

L 

-  Diffusion  length,  cm. 

e 

-  Reactivity 

i*J*.  -  Root  of  the  characteristic  equation 
■jo  -  Resonance  escape  probability 

Literature  Survey 

One  of  the  first  nuclear  scientists  to  publish  an  extensive 
derivation  of  the  "pile-kinetic"  equations  was  Hurwitz  (11).  He 
described  the  methods  used  to  solve  the  pile-kinetic  equations 
for  step  changes  in  reactivity  and  when  reactivity  was  a  slowly 
varying  function  of  time. 

Kimel,  et.  al.,  (1I4.)  presented  a  theoretical  development 
for  the  time  behavior  of  neutron  density  as  a  function  of  step 
changes  in  reactivity.  The  space  independent  kinetic  equations 
were  based  on  the  Fermi  continuous  slowing  down  model  and  used 
six  groups  of  delayed  neutrons •  Experimental  measurements  were 
made  with  the  Argonaut  Reactor  at  Argonne  National  Laboratory. 
Reactivity  worths  of  the  control  rods  of  the  Argonaut  were  meas- 
ured experimentally  using  the  theoretical  development.  Both 
positive  period  measurements  and  negative  rod  drops  measurements 
were  made  and  good  agreement  was  obtained  between  the  two  meth- 
ods, indicating  that  the  theoretical  development  could  be  used 
with  the  Argonaut  reactor. 

Smets  (17)  presented  a  solution  for  the  one  energy  group 
bare  reactor  kinetic  equations  when  the  reactivity  was  a  func- 
tion of  time.  The  time  dependence  of  the  reactivity  was  assumed 
to  be  either  linear,  exponential,  or  a  reciprocal.  Smets  com- 
bined the  basic  space  independent  differential  equations  in  the 
form 


4— •  D  + 


pfOtA, 


n(t)  = 


^,„(tVnH)     (i) 


where 


J  s  prompt  neutron  lifetime, 

Q  -   fraction  of  total  neutrons  that  are  delayed, 

A  s  fraction  of  neutrons  in  the  ith  delayed  group, 

A  s  neutron  density,  neutron/cm-', 

1  s  decay  constant  of  precursor  of  the  delayed  neutrons  of 
i 

the  ith  group, 

keff  s  effective  multiplication  constant, 
m  -  number  of  delayed  neutron  groups, 
and  D,  in  this  case  was  the  differential  operator  d/dt. 
Solutions  of  Eq,  (1)  in  the  form  of  Laplace  type  contour  inte- 
grals were  obtained.  Smets  also  found  that  the  kinetic  equa- 
tions could  be  integrated  for  any  piecewise  analytic  function, 
keff  W*  ^ich  consisted  of  straight  lines,  exponentials,  and 
inverse  functions  but  not  the  sum  of  any  two  or  more  such  func- 
tions. 

Toppel  (20)  discussed  the  errors  in  determining  the  magni- 
tude of  a  step  change  in  reactivity  by  means  of  measuring  the 
stable  reactor  period.  He  found  that  for  certain  reactivities 
the  neutron  population  would  not  be  represented  by  a  single  ex- 
ponential until  several  minutes  after  a  positive  step  was  made. 
Toppel  also  considered  a  reactor  with  a  source  and  showed  that 
a  longer  wait  time  was  required  for  the  true  exponential  neutron 
density  distribution  than  with  a  source  free  reactor.  A  typical 


point  from  Toppel's  work  showed  that  in  the  oaae  of  a  positive 
insertion  of  reactivity,   £  -  5.0  x  X0~^,   and  when  the  initial 
effective  multiplication  constant  was  1.000,  100  seconds  were 
required  for  the  reactor  to  be  within  one  per  cent  of  being  on  a 
stable  period.  In  contrast,  if  the  Initial  reactivity  was 
0,99995  and  then  the  reactivity  waa  inserted,  a  wait  time  of 
approximately  300  seconds  was  required  for  the  reactor  to  be 
within  one  per  cent  of  being  on  a  stable  period. 

Friedman  (5)  discussed  the  use  of  an  "effective"  fraction 
of  delayed  neutrons  in  place  of  the  true  fraction  of  neutrons. 
He  sts.tod  that  the  fraction  of  delayed  neutrons  should  be  in- 
creased because  the  delayed  neutrons  have  a  larger  nonleakage 
probability  than  prompt  fission  neutrons,  thus  a  delayed  neutron 
is  more  effective  in  producing  another  fission  than  is  a  prompt 
neutron.  Friedman  went  on  to  say  that  if  u238  waa  pr9gent  in 
the  reactor  fuel,  the  fraction  of  delayed  neutrons  should  be  in- 
creased due  to  the  increased  fraction  of  delayed  neutrons  from 
U2^  .  An  expression  was  given  for  the  effective  fraction  of 
delayed  neutrons  for  the  ith  precursor  group  (5). 
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.28 


(2) 


Whan 


* 


(25) 
(28) 


fraction  of  neutrons  which  are  from  the  i-th 

(25) 
precursor  due  to  fission  from  U  (28)* 


(delayed) 
(prompt  ) 


snleakage  probability  of  /d®la7®d]  neutrons 


(28) 

■NJ  (25) 
V     =  average  number  of  neutrons  per  fission  from 

(28) 
(25) 

C**  _  U26  fissions 

6    =  -§5 

U   fissions 

€.  =  fast  fission  factor  of  U238 

Ewen  and  Wick  (ij.)  extended  the  work  of  Friedman  (5)  per- 
taining to  evaluating  the  effective  fractions  of  delayed  neu- 
trons. They  gave  results  of  several  cores  and  showed  the  dif- 
ference in  core  reactivity  evaluation  by  using  their  calculated 
effective  fraction  of  delayed  neutrons.  Keepin,  et.  al.,  (13) 
measured  the  period,  relative  abundance,  and  absolute  yield  of 
delayed  neutrons  from  "fast"  fission  for  six  nuclides  including 
U  ■'"'  and  from  thermal  fission  for  throe  nuclides  including  U  -. 
The  "Godiva"  reactor  was  used  as  the  neutron  source.  Six  expo- 
nential periods  were  found  to  be  necessary  and  sufficient  for  an 
optimum  least-square  fit  of  the  experimental  data.  Keepin' s 
work  on  determination  of  delayed  neutron  parameter  is  the  latest 
presented  and  is  probably  the  most  reliable  analytical  study 
that  has  been  made  on  the  delayed  neutron  phenomenon  to  date. 
Many  persons  have  experimentally  measured  sets  of  delayed  neu- 
tron data  and  most  of  the  sets  are  slightly  different.  Keepin 
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(12)  presented  a  list  of  experimentally  measured  sets  of  delayed 
neutron  constants  and  presentod  an  excellent  description  of  the 
delayed  neutron  phenomenon. 

Skinner  and  Cohen  (l6)  showed  that  a  two  delayed  neutron 
group  model  represented,  to  a  fair  degree  of  accuracy,  the  ki- 
netic behavior  of  a  reactor  when  a  positive  insertion  of  reac- 
tivity was  considered.  They  also  demonstrated  that  at  least 
three  groups  must  be  employed  for  the  delayed  neutron  model  used 
to  describe  the  time  behavior  of  neutron  density  in  a  reactor 
when  negative  reactivities  are  inserted  into  the  core.  The 
principal  advantages  given  for  using  the  reduced  delayed  neutron 
group  models  wore:  1)  the  numbor  of  differential  equations  to 
be  solved  was  reduced,  2)  the  solutions  of  the  kinetic  equations 
required  less  computer  time,  3)  control  design  was  greatly  sin- 
plified.  The  following  is  a  list  of  three  group  constants  as 
given  by  Skinner  and  Cohen  (l6). 

Table  1.  Reduced  delayed  neutron  group  constants  for  thermal 


fission 

of 

V<H 

. 

Group 

t 

Xx  (secTM 

1 

ft/„ 

1 

0.012M). 

0.033 

2 

0.03691). 

0.3l).6 

3 

O.0632 

0.621 

Most  methods  of  solving  the  roactor  kinetic  differential 
equations  include  the  assumption  that  the  noutron  flux  can  be 
separated  into  two  variables,  space  and  time.  Solutions  of  the 
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kinetic  equations  are  difficult  to  obtain  when  the  separation 
process  is  not  performed.  Garabedian  and  Leffert  (6)  solved  the 
kinetic  equations  where  it  was  assumed  that  the  spatial  distri- 
bution contained  a  time  function.  They  found,  in  the  case  exam- 
ined, that  the  space  function  reached  an  asymptotic  shape  within 
0.005  seconds  after  a  flux  perturbation  was  caused  by  a  sudden 
reactivity  change, 

THEORETICAL  DEVELOPMENT 

Since  the  theory  used  by  this  paper  is  very  similar  to  that 
found  in  references  (8),  (11),  and  (7),  the  detailed  development 
and  solution  of  the  reactor  kinetics  equations  is  found  in 
Appendix  -  A.  The  kinetics  model  and  the  associated  assumptions 
are  given  along  with  the  resulting  equations. 

A  reactor  system  is  described  by  the  time  dependent  thermal 
diffusion  equation  and  the  precursor  rate  equations.  The  as- 
sumptions used  are  listed  as  follows: 

1)  the  reactor  is  homogeneous, 

2)  the  reactor  is  bare,  i.e.,  has  no  reflector, 

3)  the  delayed  neutron  energy  spectrum  is  the  same  as  that 
of  the  prompt  fission  neutrons. 


u 

there  is  no 

fission  from  U2^8, 

5) 

the  spatial 

distribution  of  the 

neutron  flux 

may  be 

rep- 

resented  by 

the  fundamental  mode 

i  of 

the  wave 

equatii 

m, 

■*-This  assumption  need  not  be  made  ( 

see 

Appendix 

-A). 
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6)  the  delayed  neutron  precursor  spatial  distribution  is 
proportional  to  the  spatial  distribution  of  the  neutron 
flux, 

7)  the  Fermi  age  equation  may  be  used  as  a  source  of  ther- 
mal neutrons  for  the  thermal  diffusion  equations. 

The  thermal  neutron  balance  written  around  a  volume  element 
of  the  reactor  is 

D  V%  cr,«  -  &fo4  ♦  SM  ■  ^5? .  i  %!  «* 

The  precursor  rate  equation  is 

^      -  -  *i  sM  +    qp W 

The  space  and  time  variables  of  Eq.'s  (3)  and  (1^.)  are  separated 
and  the  resulting  space  independent  kinetic  equations  ar« 


(1-P^e*?-1 


-?.tt,fa+ Wa^ 


and 

By  assuming  a  solution  for  the  set  of  seven  linear  differential 
equations  of  the  form 

nio  --  Ae*  (7) 
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\      v 


(8) 


a  characteristic  equation  may  be  obtained  as 


jlc\-  e>       m 


(9) 


or 


r- 


I  +  JLu)     i 


t — '    A' 

J >    X;q3 

+  JLu)  £ .  u>+  A^ 


(10) 


The  neutron  density  as  a  function  of  time  is  expressed  as 

T 


where  the  coefficients,  A.«s,  are 


Av   =  (I  -() 


jicv-el 


(id 


(12) 


Three  groups  of  delayed  neutrons  are  used  ag  an  approxima- 
tion for  the  usual  six  groups.  This  simplifies  the  solution  of 
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the  kinetic  differential  equations  but  they  are  of  the  same  gen- 
eral form  as  Eq.'s  (7)  and  (8). 

When  reactivity  is  a  function  of  time,  the  differential 
equations  are,  in  some  instances,  very  difficult  to  solve; 
therefore,  an  approximation  may  be  made  where  it  is  assumed  that 
the  six  groups  of  delayed  neutrons  may  be  represented  by  one 
group  of  delayed  neutrons.  McPhee  (15)  developed  a  solution  to 
the  approximated  one  group  kinetic  model  and  a  similar  model  ia 
developed  in  Appendix  -  B,  McPhee <s  equations  were  developed 
for  a  linear  time  dependent  reactivity  insertion  whereas  the 
kinetic  development  given  in  Appendix  -  B,  of  this  paper,  is 
based  on  the  assumption  of  reactivity  insertion  rate  empiricized 
to 

keII.(t)  si-  A(l-e-b*n).  (13) 

A  is  a  parameter  determined  by  the  final  value  of  the  reactiv- 
ity change  and  b  and  n  are  constants  determined  empirically* 

COMPUTED  DATA 

General 

The  solutions  of  the  differential  eouations  for  several  ki- 
netics models  were  programmed  for  the  Kansas  State  University 
IBM-6£o  computer.  Several  neutron  parameter  models  were  used. 
First,  the  usual  equations  using  six  groups  of  delayed  neutrons 
were  programmed  p.nd  two  sets  of  neutron  parameters  were  used, 
Hughes'  (10)  and  Keepin's  (13) •  Next,  an  approximation  of  the 


1$ 


six  groups  of  delayed  neutrons  by  three  groups  was  made.  The 
three  delayed  neutron  parameters  were  chosen  such  that  the  ap- 
proximation would  predict  the  same  time  dependent  neutron  flux 
after  a  step  change  in  reactivity  as  did  the  six  group  model. 
Finally,  the  kinetic  equations  were  approximated  with  ono  group 
of  delayed  neutrons  and  were  solved  for  reactivity  as  a  function 
of  time. 

The  equations  of  interest  (see  Appendix  -  A)  for  the  six 
group  analysis  with  step  inputs  of  reactivity  were  the  charac- 
teristic equation, 

o  =  J^L    +       ! 

V    JLuU+  \      ft.iiAl  / 

and  the  equation  of  neutron  flux  as  a  function  of  time  normal- 
ized to  the  critical  flux  just  before  the  step  change  in  reac- 
tivity 

An  U3M-6i>0  program  was  written  to  compute  roots  of  Eq.  (9) 
and  to  compute  the  A.'s  of  Eq.  (11).  The  computed  data  were 
used  to  analyze  the  experimental  results  from  both  the  TRIGA 


reactor  and  the  Argonaut  reactor;  therefore  two  different  prompt 

-k 
neutron  lifetimes,  2.0  x  10^  seconds  for  the  Argonaut  and 

8.0  x  10  w   seconds  for  the  TRIGA  reactor,  were  used  in  the 
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computations.  The  parameter  variation  in  each  set  of  data  was 
reactivity. 

Delayed  Neutron  Parameters 

The  two  delayed  neutron  parameter  models  used  in  the  theo- 
retical computations  were  those  as  given  by  Keepin  (13)  and  by 
Hughes  (10),  The  two  sets  of  noutron  parameters  gave  slightly 
different  sets  of  data.  The  following  tables  list  the  delayed 
neutron  constants  used  in  the  analytical  computations,  A  list 
of  several  other  3ets  of  delayed  neutron  parameters  was  made  by 
Keepin  (12). 

Table  2.  Koepin's  delayed  neutron  data  for  thermal  fission  of 
TJ235. 


Group  "  Half life  (Sec.) 


A^Soc.-1) 


*I1 


22.72 
6.22 
2.30 
0.6l 
0.236 


1.26 

0.71 

0.23 

0.0Q 

0.063 

0.025 


0.0121*4 

0.03051 

0.1114 

0.3014 

1.1363 

3.0137 


0.033  -  0.003 
0.21Q  *  0.009 
O.I96  J-  0.022 
0.395  t  0.011 
0.115  *  0.009 
O.Olj-2  t   0.008 


6  8  0.00614. 


0.034 

0.009 

0.220 

0.023 

0.282 

0.017 

0.319 

0.017 

0.112 

0.011 

0.033 

0.009 
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Table  3.  Hughes*  delayed  neutron  data  for  thermal  fission  of 
Group  j  Half life  (Sec.)   j  X   (Seo."1)         1% 

1  55.6   0.2       0.01246 

2  22.0    0.2         0.03149 

3  l).5l   0.1       0.1537 
I      1.52   0.05       0.L56 

5  0.43   0.05      1.612 

6  0.05   0.02      13.86 

(5  s  0.00755 
Reactivity  Parameter 

Computed  data  were  calculated  over  wide  ranges  of  reactiv- 
ity, from  -$0,001  to  -$10.0  for  negative  reactivities  and  from 
$0,001  to  $1.20  for  positive  reactivities  (see  PLATES  II  through 
XXI). 


EXPLANATION  OP  PLATE  II 

Flux  ratio,  (DO*)/  Q)(o)  ,  vs.  step  changes 
of  negative  reactivity,  X  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity change. 

The  prompt  neutron  lifetime,   X.  a  8.0  x  10"i> 
seconds. 

Hughes'  delayed  neutron  parameters. 
Data  computed  from  Eq.  (11). 
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EXPLANATION  OF  PLATE  III 

Flux  ratio,  (pity  (p(o)  ,  vs.  step  changes 
of  negative  reactivity,  £  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity change. 

The  prozapt  neutron  lifetime,   J_  =  8.0  x  lo"'' 
seconds. 

Hughes'  delayed  neutron  parar.eters. 
Data  computed  frorr;  Eq.  (11). 
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EXPLANATION  0?  PLATE  IV 

Flux  ratio,  §IX)  / ty(p)  ,  vs.  step  changes 
of  negative  reactivity,  ^  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity change. 

The  prompt  neutron  lifetime,  |[  ;  8.0  i  10"'' 
seconds. 

Hushes*  delayed  neutron  parameters. 
Data  computed  from  Eq.  (11). 
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EXPLANATION  OP  PLATE  V 

Negative  reactivity,    P  (dollars),  vs.  roots, 
^  (seconds"-1-),  of  the  characteristic  equation 
for  step  changes  of  reactivity. 
All  roots,  u)  ,  are  negative  in  value. 
The  A  's  indicate  the  decay  constants  of  the 
neutron  precursor  groups. 
Prompt  neutron  lifetime,   8.  s  8.0  x  10" 
second. 

Hughes'  delayed  neutron  parameters. 
Data  computed  from  Eq.  (9). 
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EXPLANATION  OF  PLATE  VII 

Positive  reactivity,   ^  (dollars),  vs. 

positive  asymptotic  period,    f-  (seconds). 

Prompt  neutron  lifetime,   8.  -  8.0  x   10"5 

seconds. 

Hughes'  delayed  neutron  parameters. 

Data  computed  from  £q.  (9)» 
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EXPLANATION  OF  PLATE  VIII 

Positive  reactivity,   ^  (dollars),  vs.  roots, 

U)  's  (seconds*^)  of  the  characteristic 
equation  for  step  changes  in  reactivity. 
All  roots  are  negative  or.cept  the  indicated 
positive  root. 

The  prompt  neutron  lifetime,   J^  3  8.0  x  lo"^ 
seconds. 

Hughes'  delayed  neutron  parameters. 
Data  computed  fron  Eq.  (9). 
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EXPIAHATION  OP  PLATE  X 

Flux  ratio,  4>(A.)/(J>(c^  ,  vs.  step  changes 
of  positive  reactivity,  ^  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity step. 

The  prompt  neutron  lifetime,   8.  =  8.0  x  10"^ 
seconds. 

Hughes'  delayed  neutron  parameters. 
Data  computed  from  Eq,  (11). 
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PLATE  X 
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EXPLANATION  OP  PLATE  XI 

Plus  patio,  (p(.-fc.)/Q>(p^  ,  vs.  atep  changes 
of  negative  reactivity,  \  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity step. 

The  prompt  neutron  lifetime,   J.;  8.0  x  10"'' 
seconds. 

Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (11). 
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EXPLANATION  OP  PLATE  XII 

Flux  ratio,  (piO/CpCo)  ,  vs.  step  changes 
of  negative  reactivity,  £  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity step. 

The  prompt  neutron  lifetime,  X  s  8.0  x  lO--3 
seconds. 

Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (11), 
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EXPLANATION  OP  PLATE  XIII 

Flux  ratio,  (\>(0/(t)(°)  •  vs.  step  changes 
of  negative  reactivity,  ^  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity step. 

The  prompt  neutron  lifetime,   j^  -  8.0  x   10 
seconds. 

Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (11). 
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PLATE  XIII 


EXPLANATION  OP  PLATE  XIV 

Negative  reactivity,    ^   (dollars),  vs.  roots, 

u) 's  (seconds  ),  of  the  characteristic 
equation  for  step  changes  of  reactivity. 
All  roots  are  negative. 

The  X  's  indicate  the  decay  constants  of  the 
neutron  precursor  groups. 

Prompt  neutron  lifetime.   Sis  8.0  x  10   seconds. 
Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (9). 
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EXPLANATION  OP  PLATE  XV 
Dimensionless  coefficients,  A  »s,  of  the  flux 

J 

equation  vs.  negative  reactivity,   if  (dollars), 

fcr  step  changes  in  reactivity. 

The  prompt  neutron  lifetime,  Hz   8.0  x  lO"-* 

seconds. 

Keepin's  delayed  neutron  data. 

Data  computed  from  Eq.  (12), 


EXPLANATION  OP  PLATE  XVI 

Positive  reactivity,   £  (dollars),  vs. 
positive  asymptotic  period,    f  (seconds), 
Prompt  neutron  lifetime,   jl  a  8.0  x  10"'' 
seconds . 

Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (9). 
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EXPLANATION  OF  PLATE  XVII 

Positive  reactivity,   (J  (dollars),  vs. 
roots,  u)'s  (seconds   ),  of  the  charac- 
teristic equation  for  step  changes  of 
reactivity. 

All  roots  aro  negative  except  the  indicated 
positive  root. 

The  prompt  neutron  lifetime,  $_  a  8.0  x  10-'1 
seconds. 

The  /A.  's  are  the  decay  constants  for  the 
neutron  precursor  groups. 
Keepin's  delayed  neutron  parameters. 
Data  computed  from  Eq.  (9)« 
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EXPLANATION  OP  PLATE  XIX 

Flux  ratio,  (}>(.t)/q>(o)  ,  vs.  step  changes 
of  positive  reactivity,  ^  (dollars),  for 
various  times,  t  (seconds),  after  the  re- 
activity stop, 

ft         m'} 
The  prompt  noutron  lifetime,   X.S  6.0  x  10 

seconds. 

Koepin'a  delayed  neutron  data. 

Data  computed  from  Eq.  (11). 
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PUIS  XIX 
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Beta  Parameter 

The  total  fraction  of  delayed  neutrons  has  been  shown  to 
vary  for  different  fissionable  materials.  Keepin  (12)  published 
a  complete  list  of  the  various  fissionable  materials  and  their 
respective  fractions  of  delayed  neutrons.  Even  with  the  same 
fissionable  material,  the  effectiveness  of  the  delayed  neutrons 
to  promote  fission  has  been  found  to  differ  from  reactor  to  re- 
actor, depending  upon  the  size  and  configuration  of  each  reactor. 
Since  effective  delayed  neutron  fractions  are  generally  dif- 
ficult to  obtain,  it  was  considered  important  to  evaluate  the 
dependence  of  the  kinetic  equations  upon  the  total  fraction  of 
delayed  neutrons.  The  ratio,   C^Ct) /($>($  ,  was  calculated  for 
various  times  after  a  step  insertion  of  reactivity  using  differ- 
ent total  fractions  of  delayed  neutrons.  The  total  fractions  of 
delayed  neutrons,  @  tSf  used  ln  the  calculations  varied  from 
0.005  to  0.008,  therefore  the  P  i3  used  cover  the  range  of  ef- 
fective fractions  used  in  most  thermal  reactors  that  utilize 
U235  as  the  primary  fuel.  The  study  included  several  negative 
and  positive  reactivities  (see  Appendix  -  F). 

The  same  computations,  as  described  above,  were  recalcu- 
lated using  a  constant  reactivity  of  -0.0748  per  cent.  A  plot 
of  the  ratio,   (j>(t)/4>(o)   »  vs»  the  total  fraction  of  delayed 
neutrons  for  two  different  times  was  made  (see  PLATE  XX).  Also 
shown  on  the  same  plot  are  the  flux  ratios  as  a  function  of  @ 
when  the  reactivity  insertion  was  -$0.10.  It  was  found  that  the 
computed  flux  ratios  for  a  particular  reactivity  ($)  insertion 


EXPLANATION  OP  PLA'EE  XX 

Flux  ratio,   (b(t)/d)(.o)  »  vs«  ttle  total 
fraction  of  delayed  neutrons,  p  < 
Flux  ratio  evaluated  at  two  time  points, 
10  seconds  and  75  seconds  after  a  step 
change  in  reactivity. 

Reactivity  was  — O.O7I4.8  per  cent, 

Reactivity  was  -$0.10,  as  defined  by 

each  (3  . 
Hughes'  delayed  neutron  parameters. 
Prompt  neutron  lifetime,  JL  s  8.0  x  10   sec. 
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was  essentially  independent  of  p  ,  whereas  the  same  was  not 
true  when  the  reactivity  insertion  was  given  in  units  of  per 
cent. 

Decay  Constant  Parameter 

The  decay  constants  of  the  neutron  precursor  groups  were 
measured  experimentally  by  Keepin,  who  quoted  the  uncertainty  of 
each  of  the  decay  constants.  In  this  paper  an  evaluation  was  to 
determine  the  effect  of  varying  each  decay  constant  within  the 
experimental  uncertainty  as  reported  by  Keepin.  A  set  of  com- 
puted data  was  made  consisting  of  the  ratio,   <ptU/(^(o}  ,  for 
various  times  after  step  insertions  of  reactivity  using  each  de- 
cay constant  modified  to  its  highest  value  of  uncertainty  (see 
Appendix  -  P ) . 

Prompt  Neutron  Lifetime  Parameter 

Since  the  prompt  neutron  lifetime  varies  considerably  from 
reactor  to  reactor,  it  is  important  to  compute  the  effectiveness 
of  the  lifetime  on  the  kinetic  equations.  For  thermal  reactors, 
lifetimes  are  known  to  vary  from  approximately  3.0  x  10"-'  sec- 
onds for  certain  swimming  pool  type  reactors  to  approximately 
1.0  x  10  *  seconds  for  large  graphite  piles. 

The  flux  ratio  equation  was  solved  for  both  positive  inser- 
tions of  reactivity  and  negative  insertions  of  reactivity  using 
various  values  of  the  prompt  neutron  lifetime.  The  range  of 
lifetimes  investigated  was  from  1  x  10"^  seconds  to  1  x  10"^ 
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seconds.  Plots  were  made  of  the  flux  ratio  as  a  function  of 
time  after  a  step  insertion  of  reactivity  for  several  reactivi- 
ties and  several  lifetimes  (see  PLATES  XXI  and  XXII).  The  plots 
indicated  that  the  lifetime  parameter  is  not  an  important 
parameter  of  the  kinetic  equations  except  when  considering  large 
step  changes  of  positive  reactivity,  ^  >  $0.30. 

Three  Groups  of  Delayed  Neutrons 

The  reactor  kinetic  equations  were  solved  using  three 
groups  of  delayed  neutrons.  The  resulting  characteristic  equa- 
tion was  programmed  and  the  roots  were  computed  by  the  IBM-650. 
It  was  demonstrated  that  the  six  delayed  neutron  group  kinetic 
equations  could  be  reasonably  approximated  with  three  groups  of 
delayed  neutrons. 

It  is  difficult  to  obtain  a  set  of  three  group  constants 
analytically,  either  experimentally  or  theoretically,  which  will 
reasonably  approximate  the  kinetic  equations  using  six  delayed 
neutron  groups.  Therefore,  a  trial  and  error  method  was  used  to 
determine  the  three  group  constants.  Flux  ratios  had  been  com- 
puted for  step  reactivity  insertions  using  the  six  group  con- 
stants. The  three  group  constants  were  then  obtained  such  that 
the  flux  ratios  as  given  by  the  three  group  constants  agreed 
with  those  as  given  by  the  six  group  constants  (see  PLATES  XXIII 
and  XXIV).  A  3et  of  neutron  constants  determined  in  this  manner 
is  given  in  Table  1).. 


EXPLANATION  OP  PLATE  XXI 

Flux  ratio,   q>(t)/(\)^  ,  vs.  time  (seconds)  after 
a  negative  atop  change  of  reactivity,  (dollars). 
Hughes'  delayed  neutron  parameters. 
Two  different  lifetimes,   5.  -  1  x  10"^  seconds. 

H;  1  x  10   seconds. 
Threo  different  reactivities,   ^  s  -¥0.10 

£  =  -$0.50 


£  s  -Ci.50 


1.0 


0  9 


0.8 


0,7 


0.6 


0.5 


0.  4 


O 


0.3 


0.2 


O.IS 


" 


PLATE  XXI 


P=  0.10 


1=  10  sec. 


60  80 

time 


100 


EXPLANATION  OP  PLATE  XXII 

Flux  ratio,       (b(-0/4>(p)  i  '3»   t3jne   (seconds)  after  a 
positive   step  change  in  reactivity,        %    (dollars). 
Hughes1   delayed  neutron  parameters. 

Two  different  lifetimes  are  shown,      |:  lx  10"''  seconds. 

(   :  lx  10"'   seconds. 
Three  differene  reactivities,      £    s  50.01 
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Table  1)..  Empirically  determined  three  delayed  neutron  group 
constants. 


Group  ^.  (Sec.*1)  ft 

1  0.79  0.0033 

2  0.091  0.003l)-9 

3  0.0117  O.OOO69 


For  each  reactivity  investigated  using  the  three  group 
analysis,  the  comparisons  with  the  six  group  analysis  gave  simi- 
lar flux  ratio  deviations  as  a  function  of  time.  At  two  seconds 
after  the  step  change  in  reactivity,  the  two  models,  the  three 
group  and  the  six  group,  gave  almost  the  same  flux  ratio.  For 
times  between  six  seconds  and  20  seconds,  the  three  group  model 
gave  slightly  higher  flux  ratios.  For  times  between  20  seconds 
and  100  seconds,  the  3ix  group  model  gave  slightly  higher  flux 
ratios.  For  times  larger  than  100  seconds,  the  three  groups 
model  gave  higher  flux  ratios.  The  difference  between  the  flux 
ratios  as  given  by  the  two  models  differed  significantly  for 
times  larger  than  150  seconds.  As  the  step  input  of  reactivity 
increased,  the  difference  in  the  flux  ratios,  as  given  by  the 
two  models,  increased.  The  largest  reactivity  used  was  -llj^ 
(see  PLATE  XXIV). 

The  six  group  code  required  approximately  100  seconds  of 
IBM-650  time  to  solve  the  characteristic  equation  and  caloulate 
flux  ratios  for  15  time  points.  The  three  group  code  required 
approximately  25  seconds  to  solve  the  characteristic  equation 
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and  calculate  the  flux  ratios  for  XS  time  points.  Thus,  if  it 
becomes  necessary  to  perform  a  kinetic  analysis  using  certain 
additional  reactivity  parameters,  the  three  group  analysis  will 
greatly  simplify  the  calculations. 

EXPERIMENTAL  MEASUREMENTS 

General 

The  experimental  data  were  obtained  with  the  TRIGA  reactor 
located  at  the  John  Jay  Hopkins  Laboratory,  San  Diego,  Cali- 
fornia and  with  the  Argonaut  reactor  located  at  Argonne  National 
Laboratory.  The  TRIGA  data  were  obtained  by  the  author  and  the 
Argonaut  data  were  forwarded  to  this  author  by  Dr.  Bill  Sturm, 
Manager  of  the  Argonaut  reactor. 

Experimental  System 

Briefly,  the  TRIGA  reactor  consisted  of  a  cylindrical  core 
containing  a  lattice  of  cylindrical  fuel  elements  surrounded  by 
a  graphite  reflector  and  located  at  the  bottom  of  a  20  foot 
aluminum  tank  (see  PLATES  I  and  XXV).  The  fuel  elements  con- 
tained fissionable  material  in  the  form  of  an  alloy  of  uranium- 
zirconium  hydride  composed  of  eight  wt.  per  cent  uranium  en- 
riched to  20  per  cent  in  U  **,  The  reactor  contained  four  con- 
trol rods.  The  control  rods  used  in  the  reactivity  measurements 
were  the  shim  rod  and  the  regulating  rod.  The  positions  of  the 
rods  are  shown  in  PLATE  XXVII.  More  descriptive  information  of 
the  TRIGA  reactor  may  be  found  in  the  references  (19),  (21),  and 


EXPLANATION  OF  PLATE  XXV 

View  of  the  core  of  the  TRIGA  reactor 
showing  the  positions  of  the  ionization 
chambers  located  Just  outside  the  core 
reflector. 
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PLATE  XXV 
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EXPLANATION  OP  PLATE  XXVII 

A  sketch  of  the  TRIGA  reactor  showing  the 
basic  core  coaponents. 
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r'LATE  XXVII 
CEMTRAL  CONTROL  ROD 
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(22). 

The  reactivity  insertion  was  made  by  changing  the  position 
of  one  of  the  two  control  rods.  The  control  rods  were  composed 
of  boron-carbide  canned  in  aluminum  tubes.  The  control  guide 
tubes  may  be  seen  in  PLATE  XXVIII.  The  rods  had  a  l5x  inch 
travel  from  their  "full  out"  position  to  "full  in".  During  the 
rod  drop  experiments  or  during  a  scram,  the  rods  fell  under  the 
force  of  gravity.  The  time  required  for  the  rods  to  fall  from 
"full  out"  to  "full  in"  was  approximately  0.35  seconds. 

The  neutron  flux  level  was  measured  with  the  use  of  two 
Westinghouse  WL-6377  compensated  ion  chambers.  The  ion  chamber 
sensitivity  was  lj.  x  lO-^"  amps/nv  and  were  electrically  compen- 
sated havinn  a  compensated  gamma  sensitivity  of  about  3  x  10 
amps/roentgen/hr.  The  compensated  ion  chambers  are  shown  in 
PLATE  I. 

The  current  from  the  compensated  ion  chambers  was  monitored 
with  two  Keithly  micro-micro  ammeters.  The  ammeters  are  shown 
in  PLATE  XXIX  Just  above  the  reactor  operator's  head.  The  sig- 
nal from  the  micro-micro  ammeters  was  recorded  by  a  Midwestern 
model  608  direct-recording,  oscillograph.   (See  PLATE  XXX). 


EXPLANATION  OF  PLATE  XXVIII 

View  of  the  THIGA  core  before  positioned 

into  the  reactor  tank. 

..total  tubes  at  the  periphery  will  hold 

the  ionisation  chambers. 

The  perforated  tubes  are  toe  control  rod 

guide  tubes  • 
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EXPrRIMEHTAL  PROCEDURE 
General 

Two  experimental  methods  were  used  to  measure  the  reactiv- 
ity worth  of  a  control  rod.  The  first  method  described  Is  by 
far  more  commonly  used  In  the  nuclear  field  for  reactivity  meas- 
urement, and  this  procedure  Is  known  as  the  positive  period 
method.  The  second  method  Is  known  as  the  rod  drop  method. 

'''he  positive  period  method  consisted  of  several  well  de- 
fined procedural  steps.  First,  the  reactor  was  brought  to  a 
critical  position  at  a  low  powor.1  The  reactor  was  maintained 
at  the  critioal  position  for  approximately  eight  minutes  to  in- 
sure that  the  criterion  for  delayed  critloallty2  was  satisfied. 
A  positive  reactivity  change  was  made  by  partially  withdrawing 
one  of  the  control  rods.  The  neutron  population  In  the  core  in- 
creased as  a  function  of  time.  One  hundred  seoor.ds  after  the 
change  In  reactivity  the  time  function  of  the  neutron  density 
was  approximately  represented  by  a  single  exponential  function. 
The  time  constant  of  the  exponential  function  was  the  reciprocal 
of  the  stable  reactor  period.  The  time  required  for  the  neutron 
density  to  increase  by  a  factor  of  "e"  was  measured  and  that 


The  TRIGA  reactor  has  a  significant  temperature  coeffi- 
cient of  reactivity  and  thus  any  temperature  change  during  reac- 
tivity measurements  would  have  surely  invalidated  the  measure- 
ments. Other  than  the  heating  effeot,  power  was  not  a  factor  in 
the  reactivity  measurements. 

Delayed  oriticality  is  defined  as  the  state  of  a  nuclear 
reactor  in  which  all  the  neutrons  resulting  from  fission,  both 
delayed  and  prompt,  are  required  to  maintain  a  chain  reaction. 
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time  was  the  stable  reactor  period  for  that  particular  reactivity 
change.  Prom  an  analytical  relationship  between  the  stable  re- 
actor period  and  the  magnitude  of  the  step  change  in  reactivity, 
the  reactivity  was  established  for  the  particular  experiment. 

Instrumentation  Calibration 

Frequent  checks  were  made  to  insure  that  the  instruments 
used  were  calibrated  and  functioning  properly.  The  linearity  of 
the  compensated  ionization  chambers  was  checked  numerous  times 
to  insure  that  the  current  output  of  the  chambers  was  propor- 
tional to  the  reactor  power  level.  The  micro-micro  ammeters 
were  checked  several  times  during  each  run  to  insure  that  the 
proper  zero  position  was  being  used.  The  Midwestern  recorder 
was  calibrated  several  times  to  insure  that  the  proper  calibra- 
tion was  maintained. 

It  was  assumed,  in  the  7RI3A.  reactor,  that  the  thermal 
power  level  was  proportional  to  the  neutron  density  at  the  posi- 
tion of  the  detector.  The  power  level  calibration  was  obtained 
by  measuring  the  average  reactor  tank  water  temperature  as  a 
function  of  time  while  the  reactor  was  operating  at  a  steady 
power.  The  temperature  measurements  were  made  at  several  power 
levels.  The  calibration  Indicated  that  the  current  output  of 
the  ionization  chamber  was  truely  proportional  to  the  power 
level  regardless  of  the  detector  position. 

During  a  particular  run,  the  micro-micro  ammeter  measured 
neutron  flux  levels  over  several  decades  of  power.  Each  decade 
of  the  micro-micro  ammeter  required  a  slightly  different  zero 
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position,  thus  when  the  micro-micro  ammeter  waa  switched  from 
one  decade  to  a  lower  value  a  new  zero  position  was  made  by  an 
adjustment  on  the  ammeter. 

The  Midwestern  recorder  operated  on  a  galvanometer  princi- 
ple, i.e.,  the  current  supplied  to  the  recorder  galvanometer  was 
proportional  to  the  current  fed  to  the  ammeter.  The  galvanometer 
reflected  ultraviolet  light  onto  a  moving  photosensitive  record- 
ing paper.  The  recorder  was  calibrated  by  imposing  a  false  sig- 
nal, of  known  value,  from  the  micro-micro  ammeter  to  the  galva- 
nometer. The  deflection  of  the  false  signal  line  on  the  recorder 
paper  was  then  measured  as  a  function  of  the  magnitude  of  the 
false  signal  (see  PLATE  XXXII). 

Positive  Period  Measurements 

The  two  control  rods  that  were  calibrated  with  the  positive 
period  method  were  the  shim  rod  and  the  regulating  rod  (see 
PLATE  XXXI).  The  safety  rod  calibration  was  furnished  by  Gen- 
eral Atomic  personnel.  The  shim  and  regulating  rod  calibrations 
were  made  using  this  safety  rod  calibration. 

The  reactor  was  brought  to  a  delayed  critical  position  with 
the  safety  rod  out  of  the  core,  the  shim  rod  completely  in  the 


core,  and  the  regulating  rod  at  a  position  necessary  for  criti- 
cality.   The  positions  of  the  control  rods  at  the  critical 
position  were  recorded  in  the  reactor  log  book.  The  reactor  was 


^The  power  level  at  the  critical  position  was  generally 
around  20  watts. 


EXPLANATION  OP  PLATE  XXXI 


View  of  the  top  of  the  TRIGA  reactor  showing 

the  control  rod  drive  motors  and  the  transient 

rod  mechanism. 

S  -  shim  rod  drive  mechanism. 

R  -  regulating  rod  drive  mechanism. 
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PLATE  XXXI 


EXPLANATION  OP  PLATE  XXXII 

Calibration  curve  for  the  Midwestern  Model  608 

direct-reoording  oscillograph. 

Current,  I  (arbitrary  units)  vs.  deflection,  0 

{am.),   of  signal  line  away  from  the  reference 

line. 

K-l  and  K-2  are  recording  channels. 
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left  in  the  critical  state  for  approximately  eight  minutes  to 
insure  that  a  delayed  critioality  had  been  obtained.  The  shim 
rod  was  then  withdrawn  from  the  core  a  short  distance  to  nake 
the  reactor  super-critical.  The  new  position  of  the  shim  rod 
was  recorded  in  the  reactor  log  book.  The  reactor  was  left  in 
the  supercritical  state  for  at  least  100  seconds  to  insure  that 
the  Increase  in  power  as  a  function  of  time  was  exponential.  A 
stop  watch  was  used  to  measure  the  time  required  for  the  reactor 
power  to  increase  by  a  factor  of  1*5*  The  time  measured  was 
then  divided  by  the  natural  logarithm  of  1.5,  0.1j.055#  to  obtain 
the  stable  reactor  period.  Several  such  measurements  were  mad* 
with  oacii  reactivity  change.  The  reactor  was  brought  back  to 
the  critical  position  at  the  original  power  level  by  nraana  of 
Inserting  the  regulating  rod  Into  the  core.  A  technique  was 
used  to  initially  overoorapensate  the  reactivity  of  the  reactor 
with  the  regulating  rod.  It  consisted  of  inserting  the  regulat- 
ing rod  into  the  core  by  an  amount  greater  than  necessary  to 
bring  the  reactor  to  a  critical  position  at  the  original  power. 
The  power  level  nas  allovred  to  undershoot  the  original  power  and 
then  the  regulating  rod  was  withdrawn  from  the  core  in  small 
steps  to  obtain  the  desired  critioality.  The  purpose  of  the  un- 
dershoot In  power  level  was  to  "balance  out"  the  effects  of  the 
delayed  neutrons  such  as  to  obtain  a  delayed  critioality  in  a 
much  shorter  time.  The  new  position  of  the  regulating  rod  was 
recorded  in  the  reactor  log  book.  The  procedure  was  repeated  by 
withdrawing  the  shim  rod  another  short  distance.  The  stepwise 
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measurements  were  continued  until  the  regulating  rod  was  almost 
completely  inserted  into  the  reactor  core.  At  that  point,  the 
safety  rod  was  slowly  inserted  into  the  core  while  the  regulat- 
ing rod  was  being  withdrawn,  to  insure  that  near  criticality 
conditions  were  being  maintained.  After  the  safety  rod  was 
fully  inserted  into  the  coro,  the  new  position  of  the  regulating 
rod  was  recorded.  The  stepwise  reactivity  measurements  were 
again  continued  until  the  shim  rod  and  the  regulating  rod  were 
completely  calibrated. 

A  theoretical  relationship  between  the  reactivity  change 
and  the  resulting  stable  period  was  used  to  evaluate  the  reac- 
tivity of  each  of  the  rod  movements.  The  relationship  (see 
PLATE  XXXIII)  was  derived  by  General  Atomic  (19).  The  reactiv- 
ity changes  were  produced  by  the  shim  rod,  thus  the  reactivity 
of  each  of  the  rod  movements  was  summed  to  obtain  the  integral 
worth  of  the  shim  rod  (see  PLATE  XXXIV).  Since  the  shim  rod  re- 
activities were  compensated  with  regulating  rod,  similar  cali- 
bration curves  were  plotted  for  the  regulating  rod  (see  PLATE 
XXXV). 

Rod  Drop  Measurements 

Rod  drop  measurements  were  made  with  two  control  rods,  the 
shim  and  the  regulating  rod.  When  available,  two  linear  power 
level  channels  were  used  for  neutron  density  measurements.  At 
other  times,  the  second  compensated  ionization  chamber  was  too 
far  away  from  the  core  to  accurately  measure  the  low  power  used 
in  the  rod  drop  measurements. 


EXPLANATION  OP  PLATE  XXXIII 

Positive  reactivity,  $      (dollars),  vs. 
asymptotic  period,   t"  (seconds). 
Analytical  data  calculated  by  General  Atomic 
(19)  using  Keepings  delayed  neutron  data. 
Total  effective  fraction  of  delayed  neutrons, 
0.0078. 
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EXPLANATION  OP  PLATE  XXXIV 

Shim  control  rod  calibration  for  the  TRIGA  reactor. 
reactivity,   ^   (dollars),  vs.  control  rod  position, 
S  (units). 

—  ■* —  Measured  by  dropping  the  shim  rod  from  various 
critical  positions  to  "full  in". 

— 9 Measured  by  the  positive  period  method. 

Xeepln's  delayed  neutron  parameters. 
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EXPLANATION  OF  PLATE  XXXY 

Regulating  control  rod  calibration  for  the 

TRIGA  REACTOR. 

Reactivity,   f  (dollars),  vs.  control  rod 

position,  S  (units). 

— •—   Measured  by  dropping  the  regulating 

rod  from  various  critical  positions 

to  "full  in". 
— © —  Measured  by  the  positive  period  method 

using  increments  of  reactivity. 
Keepin's  delayed  neutron  data. 
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The  first  step  of  the  rod  drop  measurements  was  to  bring 
the  reactor  to  a  critical  configuration  with  the  rod  to  be  cali- 
brated at  the  position  desired  and  with  the  other  rods  at  posi- 
tions necessary  for  criticality.  The  reactor  was  allowed  to 
remain  at  the  critical  positions  for  a  length  of  time  necessary 
to  insure  that  the  criterion  for  delayed  criticality  had  been 
satisfied.  The  Midwestern  recorder  was  turned  on  to  allow  for 
sufficient  warm-up  time.  A  zero  adjust  check  was  made  on  the 
Keithly  micro-micro  ammeter.  The  magnet  holding  the  control  rod 
was  then  de-energized  and  the  rod  fell  freely  into  the  core. 
Since  there  was  no  means  available  to  stop  the  rod  between  the 
critical  position  and  the  "full  in"  position,  all  of  the  rod 
drop  measurements  were  made  from  the  critical  position  to  the 
"full  in"  position. 

The  scales  of  the  Keithly  micro-micro  ammeter  ware  1  x  10 
and  3  x  10*"n.  Each  time  that  the  current  from  the  ion  chamber 
became  too  low,  the  decade  setting  on  the  ammeter  was  switched 
to  a  lower  value,  therefore  a  new  zero  adjust  was  required.  The 
zero  adjust  operation  was  performed  quickly  and  very  little  re- 
cording signal  was  lost  during  the  adjusting  period. 

A  signal  proportional  to  the  power  level  of  the  reactor  was 
recorded  by  the  Midwestern  recorder.  The  power  level  was  pro- 
portional to  the  deflection  of  the  recorder  signal  line  from  a 
reference  line.  Thus,  the  deflection  of  the  signal  line  was 
measured  and  from  the  calibration  graph  (see  PLATE  XXXII),  the 
power  level  as  a  function  of  time  was  obtained.  The  measured 
power  level  was  then  normalized  to  unity  at  the  critical  power 
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level  Just  before  the  rod  drop,  thus  the  normalized  power  level 
gave  the  ratio,  §>&)/ §)(p}    .  Prom  the  theoretical  plots  (see 
PLATES  II  and  III)  the  reactivity  of  the  rod  drop  was  determined. 
Since  the  value  of  the  ratio,   <D  W  /  &  (Pi  ,  was  measured  for 
the  various  times,  two  seconds,  four  seconds,  etc.,  it  was  pos- 
sible to  obtain  the  reactivity  at  the  various  times.  The  rod 
drops  were  analyzed  with  two  sets  of  neutron  data,  Keepin's  (13) 
and  Hughes'  (10). 

After  each  rod  drop  was  completed,  the  reactor  was  brought 
back  to  the  initial  critical  power  level,  thus  all  of  the  rod 
drop  measurements  started  with  the  same  critical  power  level,  A 
critical  configuration  was  obtained  with  the  rod  to  be  cali- 
brated in  a  new  position  and  then  the  rod  drop  procedure  was  re- 
peated. Rod  drops  were  made  from  various  critical  positions  of 
the  rod  that  was  being  calibrated,  until  the  reactivity  worth  of 
the  complete  rod  was  detsrmined.  Both  the  shim  and  the  regulat- 
ing rods  were  calibrated  in  the  same  manner  as  described  above 
(see  PLATES  XXXIV  and  XXXV). 

EXPERIMENTAL  ANALYSIS 

TRIGA  Data 

As  outlined  in  the  procedure,  two  methods  of  reactivity 
measurements  were  used  to  calibrate  two  control  rods  of  the 
TRIGA  reactor.  The  control  rods  that  wero  calibrated  were  the 
shim  rod  and  the  regulating  rod. 

When  the  positive  period  method  was  used  to  measure  the 
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reactivity  worth  of  the  shim  rod,  the  rod  was  calculated  to  have 
a  complete  worth  of  $1.87,  according  to  the  General  Atomic  (19) 
analytical  data  (see  PLATE  XXXIV).  In  these  data  General  Atomic 
used  Keepin's  (13)  delayed  neutron  data  but  instead  of  the  total 
fraction  of  delayed  neutron  as  given  by  Keepin,  they  calculated 
a  total  effective  fraction  of  delayed  neutrons.  General  Atomic 
used  as  a  total  fraction  the  value  of  0.0078  as  compared  to  that 
given  by  Keepin  of  O.OOolj.,  but  as  was  shown  in  the  analytical 
data  of  this  paper  (see  Appendix  -  P),  there  was  no  significant 
difference  between  the  stable  period  relationship  and  reactivity 
when  either  0.0061;  or  0.0078  was  used  for  the  total  fraction  of 
delayed  neutrons. 

Because  of  the  difference  between  the  two  sets  of  neutron 
constants  as  given  by  Keepin  as  compared  to  those  given  by 
Hughes,  there  was  found  to  be  a  slightly  lower  complete  reactiv- 
ity worth  of  the  control  rod  when  the  rod  was  calibrated  by  the 
positive  period  method  using  Hughes'  delayed  neutron  data  (see 
PLATE  XXXVI),  The  difference  between  the  two  delayed  neutron 
models  was  most  significantly  demonstrated  when  the  two  plots  of 
reactivity   f  ,  vs.  the  stable  period,  f  ,  were  compared 
(see  PLATE  XXXVII).  The  total  worth  of  the  shim  rod  was  calcu- 
lated to  be  $1,75  when  Hughes'  delayed  neutron  data  were  used. 

Similarly,  the  regulating  rod  was  calibrated  by  the  posi- 
tive period  method  using  Keepin's  delayed  neutron  data  (see 
PLATE  XXXV)  and  Hughes'  delayed  neutron  data  (see  PLATE  XXXVIII). 
The  regulating  rod  complete  reactivity  worth  was  calculated  to 
be  ¥1,76  and  $1,61).  using  Keepin's  and  Hughes'  data  respectively. 


EXPLANATION  OP  PLATE  XXXVI 


3him  control  rod  calibration  for  the  TRIGA 

reactor. 

Reactivity,   £  (dollars),  vs.  control  rod 

position,  S  (units). 

— ■» —  Measured  by  dropping  the  shim  rod  from 

various  critical  positions  to  "full  in", 
— © —  Measured  by  the  positive  period  method 

using  increments  of  reactivity. 
s'  delayed  neutron  data. 
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EXPLANATION  OF  PLATE  XXXVII 


Positive  reactivity,   C  (dollars),  vs.  positive 
asymptotic  period,   t  (seconds), 
K  curve:  Keepings  delayed  neutron  parameters. 
II  curve:  Hughes'  delayed  neutron  parameters. 
Neutron  lifetime,  A  s  8.0  x  10-5  seconds. 
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EXPLANATION  OP  PLATE  XXXVIII 


Regulating  control  rod  calibration  for  the 

TRIGA  reactor. 

Reactivity,   ^    (dollars),  vs.  control  rod 

positions,  S   (units). 

—  -» —  Measured  lay  dropping  tlie  regulating 

rod  from  various  critical  positions 

to  "full  in". 
9 Measured  by  the  positive  period  method 

using  increments  of  reactivity. 
Hughes'  delayed  neutron  parameters. 
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An  interesting  point  was  demonstrated  in  the  control  rod 
calibration  by  the  positive  period  method.  During  the  process 
of  positive  period  measurements,  the  gaiety  rod,  a  control  rod 
of  almost  equal  magnitude  reactivity  worth  as  the  shim  rod  or 
the  regulating  rod,  was  moved  from  its  "full  out"  position  to 
its  "full  in"  position.  Thus,  the  safety  rod  was  moved  from  one 
extreme  of  its  reactivity  worth  position  to  the  other  extreme. 
The  point  of  interest  is  that  the  position  of  the  safety  rod  had 
no  effect  on  the  calibration  of  either  the  shim  rod  or  the  regu- 
lating rod.  If  the  position  of  the  safety  rod  had  created  an- 
other parameter,  "shadowing  effect",  in  the  calibration  proce- 
dure, a  shift  in  the  reactivity  calibration  curve  would  have 
been  observed.  The  shift  would  have  been  in  the  form  of  a  fi- 
nite jump  in  the  calibration  curve. 

Forty-three  rod  drop  experimeuts  were  performed  to  measure 
the  reactivity  wortn  of  the  control  rods.  The  control  rods  used 
in  the  rod  drop  experiments  were  the  shim  rod  and  tha  regulating 
rod.  In  each  case  the  control  rod  was  dropped  from  Its  critical 
position  to  its  "full  in"  position;  thus  the  reactivity  that  was 
calculated  was  the  difference  in  reactivity  of  the  two  positions. 
It  was  assumed  that  the  rod  drop  represented  a  negative  step 
change  in  reactivity.  The  time  required  for  each  control  rod  to 
travel  from  its  "full  out"  to  its  "full  in"  position  was  approx- 
imately 0.35  seconds.  Thus,  at  times  large  compared  to  0.35 
seconds,  the  reactivity  Insertion  appeared  as  a  negative  stop 
change. 

Because  of  the  finite  reactivity  insertion,  it  was  found 
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that  for  times  less  than  eight  seconds  after  the  rod  drop  was 
initiated  the  measured  neutron  flux  was  higher  than  given  by  the 
theoretically  predicted  flux.  Since  a  step  change  in  reactivity 
was  assumed  in  the  theory,  this  observed  effect  was  expected. 
The  discrepancy  between  the  analytically  determined  flux  and  the 
experimentally  measured  flux  was  greater,  for  short  times  after 
the  rod  drop,  when  larger  reactivity  insertions  were  made  (see 
PLATE  XLII  and  PLATE  XLIII).  Since  the  reactivities  were  deter- 
mined froju  the  ratio,   (f)  Ot) /(D(P)  ,  at  various  times  after  the 
change  in  reactivity,  the  apparent  reactivity,  measured  at  times 
just  after  the  rod  drop  was  made,  rere  smaller  in  aagnitude  than 
the  true  value  of  the  reactivity  change.  In  each  case,  the 
apparent  reactivity  as  a  function  of  time  would  reach  a  fairly 
stable  value  for  timos  larger  than  ten  seconds  and  then  would 
vary  somewhat  around  the  stable  value.  There  was  found  to  be  no 
uniform  shift  la  the  reactivity  values  as  given  at  large  times 
after  the  rod  drop,  although  the  reactivities  did  tend  to  drift 
away  from  the  stable  value  of  reactivity  for  tines  larger  then 
75  seconds.  In  each  run,  the  reactivities  were  averaged  for 
times  greater  than  ten  seconds  and  less  than  100  seconds  and  the 
mean  value  was  reported  as  the  measured  reactivity  of  the  rod 
drop. 

Integral  control  rod  calibrations  were  made  for  the  shim 
rod  and  the  regulating  rod  (see  PLATES  XXXVI  and  XXXVIII).  The 
total  worth  of  the  shim  rod  was  calculated  to  be  &1«97  and  $1. 84 
using  Keepin's  and  Hughes'  delayed  neutron  data  respectively. 
The  rod  drop  method  gave  reactivity  worths  which  wore  larger  in 
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magnitude  for  large  reactivities  than  did  the  positive  period 
method  (see  PLATE  XXXVIII).  The  total  worth  of  the  shim  rod  was 
found  to  be  5.1  per  cent  higher  when  measured  with  the  rod  drop 
method  than  when  measured  with  the  positive  period  method  using 
either  Keepings  or  Hughes'  delayed  neutron  data.  (See  PLATES 
XXXIV  and  XXXVI).  This  discrepancy  is  attributed  primarily  to 
the  fact  that  the  total  reactivity  worth  of  the  control  rod  was 
measured  in  one  reactivity  step  when  the  rod  drop  method  was 
used  whereas  when  the  positive  period  method  was  used,  the  total 
reactivity  worth  was  determined  by  12  incremental  steps  of  reac- 
tivity. 

Argonaut  Data 

Through  the  courtesy  of  the  Argonaut  reactor  personnel  at 
the  International  School  of  Nuclear  Science  and  Engineering, 
Argonne  National  Laboratory,  experimental  data  were  obtained 
from  the  Argonaut  reactor.  The  data  consisted  of  several  rod 
drop  experiments  with  both  a  slab  and  an  annular  core  loading. 

An  analysis  of  the  data  gave  similarly  behaved  results  as 
that  given  by  the  TRIGA.  When  the  rod  drop  experiments  were 
analyzed  using  Keepin's  data  the  evaluated  reactivity  was  always 
higher  than  that  given  by  Hughes'  data.  The  same  type  of  dis- 
agreement was  found  between  reactivity  as  measured  by  the  posi- 
tive period  method  and  reactivity  as  measured  by  the  rod  drop 
method.  At  small  reactivities,  the  two  methods  gave  substanti- 
ally the  same  results  but  at  larger  reactivities,  the  rod  drop 
method  gave  higher  values.  The  fine  rod  of  the  Argonaut  reactor 
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was  calibrated  both  by  dropping  the  rod  from  various  critical 
positions  to  "full  in"  and  then  by  dropping  the  rod  from  a  crit- 
ical position  with  the  fine  rod  in  its  "full  out"  positions  to 
various  lower  positions.  The  two  methods  of  rod  drops  gave 
approximately  the  same  results  (see  PLATES  XLI  and  XLII). 

Reactivity  as  a  Function  of  Time 

General  Atomic  (19)  measured  the  time  required  for  a  typi- 
cal TRIGA  control  rod  to  fall  from  its  "full  out"  position  to 
its  "full  in"  position.  The  time  required  was  approximately 
0.35  seconds.  The  fall  time  was  certainly  finite  when  consider- 
ing times  less  than  ten  seconds  after  initiating  the  rod  drop. 
It  was  considered  important  to  study  the  problem  of  predicting 
the  reactivity  insertion  rate  and  the  resulting  time  dependence 
of  the  neutron  flux. 

Since  the  TRIGA  core  was  located  in  a  water  bath,  the  con- 
trol rods  fell  in  a  fairly  resistive  medium  during  rod  drop 
measurements  (see  Experimental  System).  It  was  assumed  that  the 
control  rod  motion  could  be  treated  as  a  falling  body  in  a  vis- 
cous medium  where  the  resistive  force  is  proportional  to  veloc- 
ity (3).  Therefore,  the  force,  P,  acting  on  the  control  rod, 
was  assumed  to  be  of  the  form 


Where  r  is  the  proportionality  constant,  m  is  the  mass  of 
the  control  rod,  g  is  the  acceleration  due  to  gravity,  and  4? 


EXPLANATION  OP  PLATE  XXXIX 

Calibration  curve  for  the  fine  control  rod  of  the 
Argonaut  reactor. 

Reactivity,   f  (dollars),  vs.  control  rod  position,  3. 
Upper  set  of  curves:  Keepin's  delayed  neutron  data. 
— 0 —  Measured  by  dropping  rod  from  various  critical 

positions  to  full  in. 
H    Measured  by  dropping  rod  from  a  critical 

position  at  100  to  various  lower  positions. 
—  •»--  Measured  by  positive  period  method. 
Lower  set  of  curves:  Hughes'  delayed  neutron  data. 
— A —  Measured  by  dropping  rod  from  various  critical 

positions  to  full  in. 
— «. —  Measured  by  positive  period  method. 
The  prompt  neutron  lifetime,    s  2.0  x  10~^  seconds. 
Slab  core  loading. 
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EXPLANATION  OF  PLATE  XL 

Shim  control  rod  calibration  for  the  Argonaut 

reactor,  annular  loading. 

Reactivity,    ^  (dollars),  vs.  control  rod 

position,  S  (units), 

— © Measured  by  the  rod  drop  method  using 

Keepin's  delayed  neutron  parameters. 
Q Measured  by  the  rod  drop  method  using 

Hughes'  delayed  neutron  parameters. 
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is  the  velocity  of  the  control  rod,  positive  upward.  3q.  (lij.) 
was  solved  for  the  rod  displacements,  y  ,  as  a  function  of 

time,  i.e., 

r  m  -  mgt/p  +  f[|f(0)  +  fEJd  -  •"«*)  .  (15) 

Since  the  control  rods  were  dropped  from  rest,  f?(0)  was  equal 
to  zero.  All  of  the  TRIGA  control  rods  contained  approximately 
the  sane  amount  of  mass,  therefore,  the  ratio,  m/r,  was  con- 
sidered to  be  a  constant,  c  .  The  rod  position  as  a  function 
of  time  could  then  be  written  aa 

y  s  cgt  -  c2g(l  -  e"ct)  .  (16) 

The  boundary  conditions  were 
7(0)  a   0 
and  y(0.35  Sec.)  s  -15.25  inches. 

Eq.  (l6)  was  solved  for  the  constant,  o  -  1,0298  seconds,  by 
trial  and  error. 

The  rod  position  was  then  plotted  as  a  function  of  time 
(see  PLATE  XLI).  Prom  the  rod  calibration  curve  and  the  rod 
position  curve,  a  reactivity  as  a  function  of  time  relationship 
was  obtained  (see  PLATE  XLI).  It  was  desirable  to  obtain  an 
analytic  function  for  reactivity  as  a  function  of  time  because 
most  methods  for  solving  the  kinetic  equations,  where  keff(t)  is 
a  function  of  time,  demand  well  behaved  time  functions  of  reac- 
tivity (1).  The  form 

keff (t)  a  1  -  A(l  -  e-bt  )  ,         tfc  f  (17) 


EXPLANATION  OF  PLATE  XLI 

S  curve:  Pal?,  of  a  control  rod  in  the  TRIGA 
reactor,  Eq.  (l6). 

The  position  of  a  control  rod,  S  (inches),  vs. 
tine  (seconds). 

The  zero  position  of  the  control  rod  is  the 
"full  out"  position. 

Reactivity  curve:  The  reactivity  inser- 
tion due  to  a  falling  control  rod  in  the  TRIGA 
reactor  as  predicted  by 

P=  i  a . *-btn)  . 

Reactivity,  f  (dollars)  vs.  time  (seconds). 

— g_-  where  n  is  equal  to  2.0. 

@ whore  n  is  equal  to  3«0» 

■ where  n  is  equal  to  1|.0. 
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la  used  here.  Note  that  kgfffO)  is  equal  to  unity  and  koff(t) 
Is  approximately  equal  to  (1.0  -  A)  at  large  times. 

In  order  to  solve  for  b  In  Eq.  (17),  an  arbitrary  bound- 
ary condition  was  imposed  upon  Eq.  (17)  such  that  the  time  de- 
pendent term  of  k  fAt)   bo  0.95  of  its  final  value  at  time  equal 
to  0.35  seoonds,  i.e.,  keff(0.35)  s  1.00  -  0.95A.  Corresponding 
to  various  assumed  values  of  n  the  following  values  for  b 
were  calculated: 


n  ■  2 

5   b  S 

A.$ 

n  =  3 

;  °  = 

69.9 

ft«* 

;  b  - 

300 

Since 

reactivity 

is 

e  = 

keff  " 

1  _  kex  _ 

"  keff~ 

* 

and 

?m  s 

kex 

13 

(18) 


then 

£($)  -i(i  -e-btn)  .  (19) 

P 

Reactivity,  f  ,  from  Eq.  (19)  was  plotted  as  a  function  of  time 
(see  PLATE  XLI)  for  the  various  values  of  n.  It  was  found  that 
n  equal  to  lj..0  gave  the  best  empirical  fit  to  the  reactivity 
curve. 

keff (t)  from  Eq.  (17)  was  used  in  the  reactor  kinetic  equa- 
tion (see  Appendix  -  B) .  Several  computations  were  made  and  the 
resulting  predictions  of  neutron  flux  as  a  function  of  time 
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after  a  change  in  reactivity  were  compared  to  experimentally 
measured  value3.  The  kinetic  equations  used  for  the  calcula- 
tions cere  of  the  form 


nW-nW 


a,  + 


fP-0 


Qa+C&-Me 


wr 


-i 


A<& 


><k' 


•  cu+  t^-Oe" 


Ur 


Q2e 


bT 


+  (M 

(20) 


whore  the  decay  constant,  ^  ,  represents  the  single  delayed 
neutron  group.  It  was  found  that  the  value  of  n  equal  to  lj.,0 
gave  the  best  comparison  with  the  experimental  data.  Even  then, 
although  the  predicted  flux  levels  from  the  approximated  reac- 
tivity equation  gave  better  comparisons  for  very  short  times 
after  the  initial  rod  drop  than  those  using  the  six  group  step 
input  analysis,  the  results  were  in  general  quite  poor.  It  was 
then  decided  to  let  both  b  and  ^  be  empirical  parameters. 
The  best  fit  to  experimental  measurements  was  obtained  when  (see 
PLATE  XLII) 

b  s  10.0  Sec. 
/(    -  0.20  Sec."1 

These  parameters  gave  good  results  over  a  fairly  wide  range 
of  reactivities.  The  largest  reactivity  insertion  used  was 
-$ljj.2  (see  PLATE  XLII).  The  finite  reactivity  insertion  method 
gave  decisively  better  predictions  of  the  time  behavior  of  the 
neutron  flux  for  small  times,  t  ^  four  seconds.  At  larger 
times  the  model  gave  increasingly  poorer  results  due  to  the 
assumption  of  one  group  of  delayed  neutrons. 
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CONCLUSIONS 

The  reactor  kinetic  theory  used  in  this  paper  is  quite 
satisfactory  for  use  with  the  TRIGA  reactor,  at  least  for  reac- 
tor operation  at  low  power  levels.  The  conclusion  above  is 
based  on  the  agreement  between  reactivity  determination  by  the 
method  of  positive  period  measurements  and  by  the  method  of  rod 
drop  measurements.  The  rod  drop  method  depended  upon  the  theo- 
retical prediction  of  the  flux  ratio  over  wide  ranges  of  time. 
Good  agreement  was  found  between  the  predicted  flux  ratios  and 
the  experimentally  measured  flux  ratios  for  Keopin's  and  Hughes' 
delayed  neutron  parameters.  The  agreements  between  predicted 
and  measured  flux  ratios  were  poorer  for  large  than  for  3mall 
reactivities.  Agreement  was  unacceptable  at  reactivities  greater 
than  $2.00.  For  reactivity  measurements  greater  than  .,1.00, 
there  was  observed  a  slight  scatter  of  the  experimentally  de- 
termined flux  ratios.  It  i3  not  certain  whether  the  scatter  of 
the  observed  data  was  due  to  the  inadequacy  of  the  kinetics 
model  or  to  the  failure  of  the  experimental  equipment  to  accur- 
ately measure  the  rapidly  falling  flux  levels.  A  possibility 
that  the  kinetics  model  cannot  be  extended  to  cover  large  nega- 
tive changes  in  reactivity  was  indicated  by  the  consistent  pre- 
diction of  larger  values  of  reactivity  by  the  rod  drop  method  as 
compared  with  tne  positive  period  method.  It  must  be  remembered 
that  all  of  the  positive  period  measurements  were  composed  of 
small  positive  changes  in  reactivity.  It  is  recommended  that 
when  using  the  rod  drop  method,  for  evaluating  the  reactivity 
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worth  of  a  control  rod,  the  measurements  be  made  using  moder- 
ately small  changes,  I.e.,  that  the  total  reactivity  worth  of  a 
control  rod  be  measured  by  a  series  of  incremental  rod  drops  and 
not  be  dropping  tho  control  rod  from  "full  out"  to  "full  in". 

It  was  found  that  the  two  delayed  neutron  parameter  models, 
Hughes'  and  Kespin's,  gave  similar  theoretical  predictions  of 
the  kinetic  behavior  of  a  reactor.  Keepin's  model  consistently 
gave  larger  evaluations  of  reactivity,  approximately  eight  per 
cent,  rhen  both  the  positive  period  and  the  rod  drop  methods 
were  used  to  calibrate  the  control  rods  of  the  TRIGA  reactor  and 
the  Argonaut  reactor.  The  two  delayed  neutron  models  gave  simi- 
larly shaped  control  calibrations  for  both  the  TRIGA  reactor  and 
the  Argonaut  reactor,  that  is,  when  evaluating  large  negative 
changes  in  reactivity  the  disagreement  between  the  positive 
period  and  tho  rod  drop  method  evaluations  were  similar  for  both 
reactors  and  for  both  delayed  neutron  models.  So  conclusion 
could  be  made  as  to  which  model  fit  the  experimental  measure- 
ments best,  because  both  models  wore  consistent,  within  them- 
selves, when  considering  small  changes  in  reactivity  and  both 
models  gave  similar  disagreements  for  large  reactivities.  But 
it  is  the  opinion  of  this  author  that  either  Keepin's  or  Hughes' 
model  may  be  used  in  reactivity  evaluation. 

It  was  found,  in  the  Eeta  Parameter  analysis,  that  if  reac- 
tivities are  quoted  In  units  of  dollars  there  is  no  confusion  as 
to  the  response  of  a  reactor  to  a  particular  change  in  reactiv- 
ity (&),  that  is,  regardless  of  the  total  fraction  of  delayed 
neutrons,  the  same  reactivity  ($)  Input  gave  approximately  the 
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same  predicted  flux  ratios. 

The  kinetic  equation  parameter  analysis  revealed  that  the 
only  significant  parameter  is  reactivity.  The  other  parameters 
nay  be  varied  to  any  reasonable  values  and  the  resulting  theo- 
retical change  of  the  kinetic  behavior  of  a  reactor  will  be 
slight.  The  parameter  study  confirms  a  generalized  study  by 
Henry  (9). 

The  three  group  kinetic  equations  used  in  this  paper  con- 
firm work  done  by  Skinner  and  Cohen  (l6),  in  that  the  six  group 
equations  may  be  presented  by  the  reduced  group  equations.  The 
real  usefulness  of  the  three  group  equations  probably  does  not 
lie  in  their  ability  to  predict  the  reactivity  worth  of  reactor 
control  rods,  but  rather,  to  predict  the  general  nature  of  the 
time  dependency  of  the  neutron  flux  in  a  reactor  where  a  series 
of  reactor  calculations  must  be  made  considering  many  parameter* 
in  addition  to  the  usual  parameters.  The  computer  code  devel- 
oped for  tills  paper  may  be  altered  to  consider  a  wide  class  of 
kinetic  problems.  These  computer  calculations  s.re  more  rapid 
and  easier  to  analyze  than  the  calculations  obtained  with  the 
six  group  aquations . 

When  considering  rod  drop  experiments  or  other  reactivity 
insertion  problems,  where  the  reactivity  insertion  is  rapid  but 
still  finite,  the  six  group  kinetic  equations  will  predict  the 
time  dependency  of  the  neutron  density  when  evaluated  at  times 
much  larger  in  magnitude  than  the  total  reactivity  insertion 
time.  But,  when  neutron  densities  are  evaluated  at  times  compa- 
rable to  the  total  insertion  time,  it  is  necessary  to  consider 
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reactivity  as  a  function  of  time  instead  of  to  assume  the  step 
function.  A  set  of  equations  developed  in  this  paper  approxi- 
mated the  reactivity  insertion  rate  o±   a  control  rod  drop  and 
the  resulting  neutron  flux  in  the  TRIGA  reactor.  The  predicted 
time  dependency  of  the  neutron  flux  gavo  values  which  were  in 
good  agreement  with  the  experimentally  measured  neutron  flux. 
The  approximated  reactivity  equations  gave  better  predictions  of 
the  neutron  flux  at  times  less  than  four  seconds,  after  initiat- 
ing the  reactivity  change,  than  did  the  six  group  equations 
using  a  step  input  of  reactivity.  For  times  greater  than  four 
seconds,  the  neutron  flux  was  predicted  to  a  good  degree  of 
accuracy  by  the  kinetic  equations  using  six  delayed  groups  and  a 
step  input.  The  reason  for  the  poor  results,  using  the  approxi- 
mated reactivity  equations  at  large  times  was,  no  doubt,  due  to 
the  one  delayed  neutron  group  assumption. 
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Part  -  A 

Theoretical  Development  of  the  Kinetic  Equations 

The  system  considered  was  a  bare,  homogeneous,  thermal 
reactor.  A  thermal  neutron  balance  on  a  volume  element  of  the 
reactor  core  gives  (8) 

DV^CrO-lJEnrr-u  +  Scc.t)--  ^ni£&.  L  jfesa   (20) 

The  Fermi  age  equation  is  used  as  the  source  of  neutrons  for  the 
thermal  diffusion  equation.  The  source  term  for  prompt  neutrons 
is 

Vro**  pt  X 

The  rate  equation  for  the  precursor  density  la 

KvfcjO   a  _ ^ c       4  kjst^W  (22) 

it  f. 

where  -^.-,C;  gives  the  loss  of  neutron  precursors  due  to  radio- 
active  decay  and  ~Z^\i.^^S\\)      gives  the  production  of  precur- 
sors from  the  fission  process.  Since  delayed  neutrons  are  pro- 
duced following  the  beta  decay  of  the  precursors,  the  total 
production  rate  of  all  delayed  neutrons  is   )     AC-    •  Since 
there  are  six  groups  of  delayed  neutrons,  m  is  equal  to  six  in 
the  above  summation.  Delayed  neutrons  are  born  with  energies 
between  0.2  and  0.8  Mev.,  thus  the  production  rats  of  delayed 
neutron  term  must  be  multiplied  by  the  nonleakage  probability 


132 


during  slotting  down  and  by  the  resonance  escape  probability. 


r  ^      -B  c. 

b(r.O        =     >    -f-.e.        XC(ri\ 


-tft, 

(23) 


The  source  term  for  both  the  prompt  neutrons  and  uhe  de- 
layed neutrons  is  substituted  into  the  thermal  diffusion  equa- 
tions and  the  resulting  equation  Is 

<2W 


-t- 


£T,eBW^.^L 


Tho  following,  well  known,  relations  were  substituted  into  Eq, 
(2«4.) 


(25) 


*-.  '/{.V  (26) 

[  2J+ )  was  then  divided  by   £  ,  which  gives 

\_ 
I 


The  space  and  time  variables  of  the  thermal  diffusion  and  the 


(27) 


precursor  rate  equations  were  assumed  to  be  separable,  that  is 
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$(r,U   -  Ftr>#C«  (28) 

and 

c-<x,^=  kctitt&i .  (29) 

It  is  assumed  that  the  spatial  distribution  of  the  thermal  neu- 
trons can  be  given  by  the  fundamental  mode  of  the  wave  equation, 

V2F(r)  +  B*F(_n  -  O  .  (30) 

2  2 

Since  V   is  a  spatial  operator,   B  F(r^  is  substituted  into 
Eq.  (27)  after  the  separation  of  variables,  thus 

L  gt"  ,         (31) 

The  definition  of  the  effective  multiplication  factor  and  the 
finite  prompt  neutron  lifetime  is  used  to  simplify  Eq.  (32) 

Z  .,  =■   r~?   -   effective  multiplication  factor   (33) 

e+*  I  +•  LB 


o     -      *- 

X,     ~  =        finite  prompt  neutron  lifetime         (3I4.) 

l4l?B 
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Eq.  (32)  is  divided  by   |  +  L^g2"  and  then  Eq.  (33)  and  Eq, 

( 3^4- )  is  substituted  into  Eq.  (32). 

Eq.  (3D  is  then  divided  by  the  spatial  flux  function  to  give 


+ 


£*(i4-l5ehFtr) 

1=  1 

Eq.  (35)  is  divided  by  the  neutron  velocity,  or  ,  and  the 
expression 

iM     -    n(-t) 

V 
is   substituted  into  Eq.    (35)   giving 

I  V-1   -BT 


The  relationship 


£    * 


1 


louirC  I  V  LaB*} 


5"rr       =      finite  prompt  neutron  lifetime 


is   substituted  into  Eq.    (36)    to  give 

_   cW+)  (37) 

at 


13$ 


-tftj 

— P - — - —      is  defined  as  the  space  independent 

F«r) 

constant.  £.  . 

The  precursor  rate  equation  is  then  written  in  the  separated 

variable  form  as 

The  resulting  equation,  (39).  contains  functions  of  time  alon6 
on  the  left  side  and  functions  of  space  alone  on  the  right  side. 

dH,(^    i    +  *,h^)   ._  feJU^fjcj  (39) 

Therefore   fc;   is  a  time  independent  parameter. 

The  space  independent  neutron  density  equation  is  then 
written  as 


The  space  independent  precursor  equation  is  then  obtained 
by  dividing  Eq.  tipQ:)  by  G^-Cr),  which  gives 

The  last  term  on  the  right  hand  side  of  Eq.  []£)  is  then  multi- 
plied and  divided  by  \reB   (l  +  lfB*)  to  give 
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dHi (*)__,  „.,„.  feeBT 


£..    then  defined  as — =rrr        ,  thus  (J..      Is 

similar  to  £   of  Eq.  (35)  and  likewise  time  independent.  The 
resultant  space  independent  precursor  rate  equation  is  obtained 
aa 

The  usual  definition  of  reactivity,   £  ,  is 


Thus 

(W 

Pe«    =     ■* 

i  -e  . 

(I»5) 

£■<!•   (kS)   is   substituted  in  the  space   independent  kinetic 

equa- 

tion  to  give 

Ut-ei          Z_,                       at    . 

(M>1 

The  same  substitution,  as  above,  into  Eq.  (lj.3)  results  in 

dt  JLU~e>  . 
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Eq.'s  (I4-6 >  and  (lj.7)  constitute  a  set  of  seven  linear  differen- 
tial equations,  therefore  solutions  of  the  form 


U;  to  -  B.'  e1*  .  *9) 


f 

i 

may  be  assumed. 

Substitutions  of  Eq.»a  (I4.8 }  and  (I4.9)  into  Eq.'s  (!+6)  and  (!j.7) 

give 


(■°-tii*-,»)A'  +  E£iAiBi''0  (50) 


1=1 


and 

1 


-  cuj+a,)b;  -  0 


Ui-e) 


(51) 


A  solution  for  the  coefficient  A»  is  displayed  by  Cramer's  rule 
as  follows: 
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A 


o 
o 


0 
0 


*.e.      ¥t     \\     \S  A#*    V* 

(uo+x.)       o  o  o  o  o 

o     -(w+AO      0  o  o 

o  o     -(uhx3)      o  a 

o  o  o     -(w+^  o 


o 
o 
o 
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0 


o      -  (w+At,) 


M£|){  to   to.    \s    \s 


^«»- 


*A 


=  o 
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Since  A«  is  not  identically  zero  except  for  the  trivial  solu- 
tion, the  determinant  of  the  coefficients  must  be  identically 
equal  to  zero.  The  determinant  of  the  coefficients  was  simpli- 
fied by  multiplying  the  following  columns  by  the  indicated 
products: 

column  1  by    (U)+X|j(w+|0      ■'    (  W  +  Au^ 
column  2  by    (UM  T^UuJ^)  •  •  •  (^ +Au^1< /U  l~  e  ) 
column  3  by   CuH  AOUJ+AiV  ■  '  Cw+  A^  &£',  /£  (  V-  fc) 
column  fc  by   (u>V  *  ,KW+>OtW^>  O  "  '  "  CuHAu^  /£  Utf 
column  J  by  Cuj+^  ^  (uj  +  A^UJ+A iV^th )  (uJUb)  fe%j  1  (  ,_e) 

column  6  by  Cu»A«H«»  +  ^."  •  (W+A+KW+^yGs  /jLCl-€l 
column  7  by    (  U>+>,H  UH  A^  •  •    (uH^  ^^  /^(\-^) 

Each  of  the  columns,  except  column  1,  was  added  to  column  1, 
thus  the  determinant  is  cast  into  the  following  form. 
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The  product,  AIJKLMH,  is  set  equal  to  zero.  Since  the  terms, 
I  through  N,  are  product  terms  of  the  reactor  parameters, 
which  are  finite,  they  cannot  be  equal  to  zero,  hence  the  A 


Xij-0 


term  must  be  identically  aero.     The  tern.  A,  in  expanded  form  ia 


UO  -    - 


IM 
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Equation  (54  )   is  divided  by       TT(UJ+A1')        which  gives 

Ui-0 

hiM)] 

(55) 

Since 

then,   Eq.    (55)  becomes 

u 

(56) 

nn 


and 


JluJ+l       +    HuJ+l    /__,     uJ+A-, 
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./ 


£-,&,  Pj  from  Eq.  (56)  is  defined  as  the  effective  fraction  of 
delayed  neutrons,  i.e.  €.^  gives  the  relativo  effectiveness 

of  a  delayed  neutron  for  producing  a  sue seeding  fission  as  eom- 

i 

pared  to  a  prompt  neutron.   €.&j  *a  the  ratio  of  the  product 
of  the  resonance  escape  probability  and  the  nonloakage  proba- 
bility during  slowing  down  for  each  of  the  precursor  neutron 
groups  to  the  product  of  the  resonance  escape  probability  and 
the  nonleakage  probability  during  slowing  down  for  fission  spec- 
trum neutrons,  that  is 

■{2e«-   was  substituted  into  Kq.  (51)  and  the  resulting  equa- 
tion is  the  characteristic  equation  for  the  kinetic  differential 
equations, 

(p 

(59) 


Equation  (59)  is  a  seventh  degree  polynomial  equation  in   i*> 
and  thus  there  are  seven  roots  of  Eq.  (59) •  The  roots  are  the 


Ite 


•igen-values,  Jm1j<juis..,iu1     ,  of  the  characteristic  equation. 
Toe  solution  for  the  neutron  density  as  a  function  of  tine  is 
then  expressed  as 


^j-t  (60) 


and  the  precursor  density  is  expressed  as 


Witt)  =  £  ^fe!°j 


^jt  (61) 

y. 1 

The  two  coefficients,  A.  and  B^.,  are  not  independent  but 
amy  be  related  by  the  substitution  of  Eq.»s  (60)  and  (6l)  into 
the  rate  equation  giving 

y        }  C  \  -  ft\  X 

Therefore,  (62) 


<-'  A  • 


The  precursor  density  equations  nay  then  be  expressed  as 

H,(-n  '- r^   S"  -A'  -  -  <fc> 


1^3 


Y\. 


m 


"tirv\e  ,  "t 


-t' 


Pig.  1 


The  neutron  and  precursor  densities  are  evaluated  immedi- 
ately before  and  Immediately  after  a  step  change  in  reactivity 
for  the  following  hypothetical  case.  The  reactor  ia  brought  to 
a  critical  configuration  and  the  neutron  density  distribution  of 
time  is  considered  as  shown  in  Figure  1.  A  step  change  in  reac- 
tivity occurs  at  time,  X.  .  At  time,  "t  —  ^  ,  the  reactivity, 

P   ,  is  equal  to  zero  and  at  time,  i_  -t- ^.  ,  the  reactivity  is 
equal  to  a  constant. 

Before  the  change  in  reactivity,  the  neutron  density  is 
expressed  as 


n»  •  A.,e      +  /\.,e-      ±   •  .  •  +  Ao,e 


(«0 


When  the  reactivity  is  zero,  one  of  the  roots,  UJ0  ,   is  zero. 
The  fact  that  one  of  the  roots  is  exactly  equal  to  zero  may  be 
seen  by  writing  the  characteristic  equation  for  zero  reactivity. 


m 


SLud  \ 


e-  O  -  J^L     +  _!_  V  ^^  ■         (65) 


\-  t 


The  term,  jLiO-t-  I   ,  may  be  cancelled  from  Eq.  (65)  and  the  re- 
sulting equation  is  of  sixth  degree  in   U)  .  The  remaining 
roots,  all  negative,  have  been  called  the  "natural  vibration 
rates"  of  the  pile  by  Zoodak  (18).  Since  the  reactor  is  assumed 
to  be  critical  in  the  delayed  sense,  the  terms  of  Eq.  (614.), 
Ao,S   .  A0.e    -  ■  ■  AoC  ;  truly  vanish  leaving  only  the 

*     '     9  1  )   7 

first  term,  ^  .  The  neutron  density  is  then  a  single  con- 
stant at  time,  -^'_  £  . 

nU'-£)  -  Ac,  =  n.    .  (66) 

After  the  change  in  reactivity,  the  neutron  density  is  expressed 


(67) 


The  precursor  density  before  the  change  in  reactivity  is  evalu- 
ated from  Eq.  (62)  and  Eq.  (63)  to  be 

Hi  U '-£)  -  L^  :  ftnCV-Ofe', .  Sf/n^      1W 

*<  *,  A;    . 

Since  n(-t)  must  be  continuous  at  4l  ,  the  limit  as  £, »•  0  of 

Eq.  (6?)  is 

n(t)  =  \     A^e  r     =  n.  (69) 


i4S 


Since  "v  is  an  arbitrary  time,  a  more  convenient  time  was  chosen 

•  i 
for  the  step  change  in  reactivity  to  occur,  i.e.,  1  =  0  ,  then 

Eq«  (70)  may  be  written  as 
7 

n(o)  =   \     Ay  =  n.  (71) 

H 

The  equilibrium  precursor  density  is  then 

14.  (0)     =  ^i6'    Y\{o)  (72) 

A, 

The  Laplace  transform  of  the  spatially  independent  kinetic 
equation,  Eq.  (lj.0),  is 

L.mn.)-    Moi)     *     —    fil?\f\U04  V,>,FW)      (73) 

\  =  I 

Also,  the  Laplace  transform  of  the  precursor  equation  is 

Eq.  (7i)-)  is  rewritten  as  St  I  I  -  €  ) 


lit* 


H;(A-)  from  Eq.  (75)  is  substituted  into  Eq.  (67)  and  the 
result  is 


fontofc; 


6-MM-nCo 


Since 


P*T> 


and 


1=1 


n  (0+)  -  n. 


then  Eq.  (76)  may  be  written  as 


ni/u 


*u-v 


t«\ 


Eq.  (77)  may  be  solved  for  Y\  (.a,), 


hCM 


&. 


'MiAiii 


The  characteristic  equation  may  be  rewritten  as 


(77) 


(78) 
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(79) 
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The  denominator  of  Eq.  (78)  may  be  recognized  as  exactly  the 
same  form  as  Eq.  (79)  and  since  the  roots,  u)  ,  satisfy  the 
characteristic  equation  the  denominator  of  Eq.  (78)  Is  zero. 
Therefore,  an  inverse  transform  of  Eq.  (78)  is  of  the  form 


**'  (80) 

vO„\. 


yv^e* 


A  method,  as  given  by  Wylie  (33),  predicts  the  coefficients, 
A^lu>j)  A-  ,  by  using  the  following  form  of  the  inverse  transform 


frH-i^a*"*4 


■1 

where  CL  (a,)  is  given  by  the  equation 
d 


t 


W) 


U--m; 

a 


%%j  A;  gj  (82) 


Using  Eq.  (82),  the  coefficients  are  determined  as 


Ajiu)j)=  nco)  (i-e) 
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Therefore,  the  neutron  density  equation  is 


(i-e) 


„      t1^ 


JUi-o 


which  is  called  the  neutron  density  retio  or  flux  ratio. 

Part  -  B 
Reactivity  as  a  Function  of  Time 

When  the  coefficient,  keff,  in  Eq.'s  (lj.0)  and  (1+3),  remains 
a  constant  or  is  a  step  function  of  time,  the  kinetic  differen- 
tial equations  are  easily  solved.  When  keff  has  a  time  depend- 
ence, the  kinetic  differential  equations  become  difficult  to 
solve  exactly  (1).  The  degree  of  difficulty  depends  upon  the 
time  function  which  describes  keff(t)  and  on  the  number  of  de- 
layed neutron  groups  used  in  the  model. 

Glasstone  (7)  presents  an  approximate  solution  of  the 
spatially  independent  kinetic  equations  where  keff  is  a  linear 
function  of  time.  One  group  of  delayed  neutrons  was  assumed  in 
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Glasstone«3  solution, 

A  solution  to  the  spatially  independent  kinetic  equations 
is  given  in  this  paper  in  which  keff(t)  is  given  by  Eq.  (ll).). 
Eq.  (1I4.)  is  an  approximation  to  the  reactivity  insertion  rate 
for  a  rod  drop  in  the  TRIGA  reactor.  For  simplicity  of  solu- 
tion, it  is  assumed  that  there  is  only  one  group  of  delayed  neu- 
trons. The  usual  assumption  connected  with  the  simplified 
kinetics  model  are  made  (see  Appendix  -  A).  Also,  in  accordance 
with  an  assumption  made  by  McPhee  (15),  in  a  similar  develop- 
ment, the  prompt  neutron  lifetime  was  assumed  to  be  approxi- 
mately zero. 

The  spatially  independent  kinetic  equations  are  the  same  as 
Eq.  (lj.0)  and  Eq.  (43 )#  except  that  here  there  is  assumed  to  be 
only  one  group  of  delayed  neutrons.  Thus  Sq.'s  (1).0)  and  (2^3 ) 
become 


MiWu-i 


'-f1*  MHai 


(85) 


d_M)  .  _AHto  +  fc.Hfti"H>ij' 

2\i(i~)  is  substituted  from  Eq.  (86)  into  Eq.  (85)  to  obtain 
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(86) 
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It  was  also  assumed  that  the  effective  fraction  of  delayed 

i 
neutron  constant,  £  £,  ,  is  equal  to  unity.  Thus,  the  two 

terms  cancel  in  Eq,  (87). 
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Eq.  (86)  is  then  differentiated  with  respect  to  time  to  give 


d2n(l)       UmK)-  01    dn(.-t)  .   nfo  dL<tfr) 


di 
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Differentiating  Eq.    (86)   gives 
Si 
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Eq.  (88)  is  then  solved  for 


[w*>- 


dH(t) 
~dF* 


J  JL£    di 

QHU.)  Is  substituted  from  Eq,(91)  into  Eq.  (90)  to  give 


(91) 


dt* 


rke«(-ti  -f 


JIG     dt   g     9 
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Then,  froc  Eq.»a   (89)  and  (92), 


is  found  to  be 


d2^)  .    (le^l-U-l)  I  •  dn(fc)   4-  n^d^fe). 


dt 


dT 


5-       d-t 


*  X^WfekJ  )n(t)-^^ 
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Eq.  (93)  is  simplified  into 

dfn(t)   _  d  nft)  .  1 

dt      a 


dV 


dt 


feef+C-t)-\-AL-Wt+flt^ 
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dt  dt   .. 


Following  McPhee,  it  is  assumed  that  the  products, 


dtZ 


7\SL  are  approximately  zero,  therefore  Eq.  (9lj.)  becomes 
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Eq.    (95)  becomes 
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Terms  in  n(t)  are  collected  to  give 

-fat* 
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E<J»  (99)  is  integrated  from  zero  time  to  time  t. 
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Where      Qt  is 


(101) 


The  integration  is  performed  directly  with  one  of  the  terms  of 
Eq.  (100)  to  give 


»-iK-4bt'*)ew'  dt' 
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cu  +  l^-Oe 
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The  last  two  terms  of  Eq.  (100)  are  integrated  numerically  with 
the  IBM-650  using  Simpson's  Rule.  The  neutron  flux  is  then 
represented  by 


n(t)  =  nCo) 


ciz.  +  ie-n 
&z  +  te-ne 
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Computer  Program  for  the  Klnetlo  Equations  Kith  Six 
Groups  of  Delayed  Reutrons 

General.     This  program  solves  the  reactor  kinetic  equations  for 
a  bare  homogeneous  thermal  reactor  assuming  that  there  are  six 
delayed  neutron  groups   (Appendix,   Part  -  A).     The  code  was 
written  In  SOAP-II  for  the  IBM-650  computer.     The  output  data 
from  the  code  produces  the  following  quantities: 

1)  Coefficients  of  the  characteristic  polynomial. 

2)  Roots  of  the  characteristic  equation. 

3)  Coefficients  of  the  neutron  flux  ratio  equation, 

k)     The  ratio,    ty&)/cjpG>)  ,  for  various  times  after  a  step 
insertion  of  reactivity. 
The  Coefficients  of  the  Characteristic  Equations.     The  expansion 
of  the1  characteristic  equation  may  be  written  as 


L  -I         uj+Ai       uj  +  1-l  uu+Au 

By  multiplying  both  sides  of  Eq.   (97)  by        ~|J".    (uj-t-V") 
and  collecting  terms  in  powers  of     U)  ,  the  following  type 
equation  results 
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Roots  of  the  Characteristic  Equation.  Once  the  reactor 
parameters,  jl,   p,    J\.&,  and  b-'n,   are  known,  the  numerical 
coefficients  of  the  characteristic  equation  may  be  calculated 
and  thus  the  particular  characteristic  equation  is  determined. 
Roots  of  such  an  equation  may  be  determined  in  many  different 
fashions .  Since  all  of  the  roots  of  the  equation  are  real  and 
the  function,  4(uj)»  is  differentiable,  the  Newton-Raphson 
method  may  be  used.  The  Newton-Raphson  method  is  an  iterative 
procedure  and  convergence  to  a  root  of  the  polynomial  equation 
is  fast  if  the  equation  is  well  behaved. 

Educated  guesses  for  the  roots  are  aided  by  the  fact  that 
the  characteristic  equation  has  six  poles,  located  at  the  nega- 
tive value  of  each  of  the  decay  constants  of  the  precursor 
groups.  Test  for  proper  convergence  is  simple  since  the  roots 
must  lie  in  well  defined  positions.  One  root  must  be  algebrai- 
cally larger  than  the  negative  value  of  the  smallest  decay  con- 
stant and  one  root  must  be  algebraically 'smaller  than  the  nega- 
tive value  of  the  largest  decay  constant.  Each  of  the  other 
roots  lie  between  respective  poles. 

For  the  roots  to  lie  between  two  poles,  one  of  the  pole3  is 
tried  for  a  guess  of  the  root.  If  proper  convergence  is  not 
obtained,  the  other  pole  is  tried  as  a  guess  for  the  root.  If 
proper  convergence  is  still  not  obtained,  the  midpoint  between 
the  two  poles  is  tried  as  a  guess.  Prom  experience  with  the 
program,  the  three  guesses  are  sufficient  to  extract  the  desired 
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root  when  tho  root  lies  batneea  too  poles.  Since  the  largest 
negative  root  io  strongly  dependent  upon  the  magnitude  of  the 
step  change  in  reactivity,  a  series  of  guesses  are  offered. 
First,  tho  largest  negative  pole  is  tried  as  a  guess  and  if  con- 
vergence is  not  obtained  to  the  proper  value,  a  multiple  of  the 
pole  magnitude  is  tried.  The  largest  algebrr.ic  root  is  always 
algebraically  larger  than  the  negative  value  of  the  smallest 
decay  constant.  If  the  reactivity  insertion  is  negative,  this 
root  will  be  negative  and  if  the  reactivity  insertion  is  posi- 
tive, this  root  will  be  positive. 

If  the  reactivity  insertion  is  negative  several  guesses  are 
tried,  i.e.,  sero,  the  largest  pole,  and  the  midpoint  between 
zero  and  the  largest  polo.  If  the  reactivity  is  positive  an 
approximation,  as  given  by  Glasstone  (7),  is  used  as  a  guess  for 
the  root. 

1,1    _     A  (XV  6 

-  —D —  (105) 

Where  \MJi   is  the  average  decay  constant  and  is  defined  as  the 
inverse  of  average  mean  life  of  the  delayed  neutrons,  i.e., 

.     _   J 

Ao.-/  —   t.    |  (106) 

\  -- 1      A  \ 

If  convergence  to  the  desired  root  is  not  obtained,  a  multiple 
of  the  approximation  is  used. 

Coefficients  of  the  Flux  Equations.  The  coefficients,  Ai, 
of  the  flux  equation 
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A; 


ui-fc 


e. 
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*  (107) 

r-1 

as  defined  by  Eq.  (9^)  are  computed  In  a  straight  forward  man- 
ner. Since  this  is  very  simple,  a  discussion  of  the  procedure 
is  not  given. 

Flux  Ratios.  Since  all  of  the   U)  'a  and  A,»s  have  been  com- 
puted, it  is  a  simple  manner  to  calculate  the  flux  ratios. 
Program  Operation.  If  one  wishes  to  compute  continuously  a  set 
of  varying  parameters,  such  as  reactivity,  the  following  two 
cards  are  added  to  the  program  just  before  the  parameter  data. 

Symbolic  Instruction   Op   Data   Address 

HZC  CON65  C0H80      1683      k°       I633    153^ 
C0N80      NOP  0000  8000       151)_3       00   0000    8000 

The  position,  l53kt   wlll  tb-en  °e  removed  from  the  table  of 
availability.  The  various  parameter  changes  may  be  added  to  the 
end  of  the  data  and  a  transfer  card  positioned  between  each  data 
card. 

The  computer  will  tend  to  overflow  when  larger  positive 
reactivities  are  used,  ^  >  $1.00.  The  overflow  is  due  to  the 
unruly  nature  of  the  characteristic  polynomial  at  larger  reac- 
tivities and  no  special  programming  precaution  has  been  used  to 
prevent  the  overflow  in  this  case .  In  order  to  prevent  large 
exponential  values  from  causing  the  computer  to  overflow,  all 
values  "e"  to  the  power  2  135.0  and  larger  are  omitted.  At 
larger  reactivities,  ^  7  §0.6o,  the  computer  will  give  errone- 
ous results  for  large  times  after  the  step  change  in  reactivity. 
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to> 


27 


This  is  no  real  handicap  since  ~i 2.  |Q  has  no  practical 

value  for  the  usual  reactor  operation. 

Input  Data.  The  nuclear  parameters  are  placed  in  the  code  in 

the  following  positions. 


(smallest  decay 
constant) 
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Aw 
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RHO 

? 

0078 

L 

Prompt  Nei 

(the  fraction 
associated  with 


Au> 


Output  DATA.  The  output  data,   U)V  *  H's,  A1«s,  and  1-—  'flare 

given  in  tha  following  ordar: 

First  line  of  data  -  all  of  the  decay  constants  and  the  average 

decay  constant. 

Second  line  -  All  of  the  delayed  neutron  fractions  and  the  total 

fraction. 
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Third  line  -  Reactivity  and  the  prompt  neutron  lifetime. 

Fourth  line  -  The  coefficients  of  the  characteristic  polynomial. 

Fifth  line  -  Hoots  of  the  characteristic  equation. 

Sixth  line  -  Coefficients  of  the  flux  equation. 

Lines  seven  through  22  -  Flux  ratios  at  the  various  times. 

After  the  data  cards  have  been  punched  and  positioned  in 
the  program  deck,  the  following  console  settings  are  used. 
Storage  Entry  Switches      70  1951  1999     Plus 


Programmed  Stop 

STOP 

Half  Cycle 

RUN 

Control 

ROT 

Address   Selection 

Anything 

Display 

Anything 

Overflow 

SENSE 

ERROR 

STOP 

EXPLAHATION  OP  PLATE  XLIV 

Plow  diagram  of  the  computer  program 
for  the  solution  of  the  reactor  kinetic 
equations  using  six  groups  of  delayed 

neutrons . 
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Part  -  D 

Three  Delayed  Neutron  Group  Program 

General .  The  basic  principles  of  the  computer  program  used  to 
solve  the  three  group  kinetic  equations  were  the  aaae  as  those 
used  in  solving  the  equations  for  six  delayed  neutron  groups. 
The  program  was  simplified  considerably  since  the  coefficients 
of  the  characteristic  polynomial  were  greatly  simplified.  Since 
the  roots  of  the  characteristic  equation  were  farther  apart,  in 
magnitude,  additional  trials  for  the  roots  between  the  poles  of 
the  equation  were  required  in  order  that  the  Newton-Raphson 
method  would  converge  to  the  proper  root. 
Input  Data.  The  required  input  data  is: 
Position        Parameter 

1387  /.?   (smallest  decay  constant) 

1388  ~A-l 

1389       x, 

1377  6   (Associated  with  the 

3   smallest  decay  constant) 

1378  ^ 

1379  p, 

Hfcl         Reactivity,  £ 

0028         Prompt  neutron  lifetime,  jl 

0236  Average  decay  constant,  /^ 

0020  Total  fraction  of  delayed 

neutrons ,  Q 

Output  DATA,  The  output  data  appears  in  the  same  order  as  in 

the  3lx  group  program. 
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Part  -  E 


Reactivity  as  a  Function  of  Time 


The  equation  that  required  a  numerical  integration  was 
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The  two  integrals  of  Eq.  (108)  were  solved  numerically  by  use  of 
Simpson's  Rule, 


(109) 


where  the  interval  from  0 
Thu«   . 

-t-o 


t  was  divided  into  2m  increments. 


ln  = 


zo 
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(110) 


The  two  integrals  of  Eq.  (108)  are  almost  alike,  thus  the 
same  part  of  the  program  was  used  to  evaluate  the  two  integrals. 
The  two  integrals  were  identified  by  index  register  A.  Since 
the  integrals  were  evaluated  at  various  times  after  the  reactiv- 
ity change,  time  was  indexed  by  register  C.  The  integrals  were 
evaluated  between  time  (C)  and  time  (C  -  1)  and  then  added  to 
the  value  of  the  integral  from  time  (0)  to  time  (C ) . 
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A  continuous  looping  system  was  devised  such  that  the  com- 
puter kept  track  of  the  point  being  calculated  to  determine 
whether  it  was  an  odd  or  an  even  time  point.  If  the  point  was 
odd,  the  value  of  the  function  at  that  point  was  multiplied  by- 
four  and  if  the  point  was  even,  the  function  evaluated  at  that 
point  was  multiplied  by  two,  as  indicated  by  Eq.  (109).  After 
the  function  was  evaluated  at  a  particular  time,  the  next  time 
point  was  obtained  by  adding     to  the  last  time  point.  The 
two  end  points  y0  and  ypn,  were  determined  by  testing  the  time  to 
see  if  It  corresponded  to  the  end  point  times. 

Since  the  exponential  function  grew  very  large  or  very 
small  as  time  increased,  a  test  was  used  to  see  whether  the  ex- 
ponential was  larger  or  smaller  than  135.0.  An  overflow  by  the 
computer  would  have  destroyed  the  results  of  the  computation. 

The  following  symbolic  terms  were  used  In  the  programt 


SUM 


I  c 


PART     -        ^0     +-  4- ly,   +2i|2+4Hj+  2.1*4.-*-  -   •  ■■»■  l/L 

C_=0     I  c 
NOPT  -  Neutron  density  as  a  function  of  time. 

T     -  Time. 

A2    -    I  _  0  +  £ 

A     -  Magnitude  of  the  final  value  of  kex« 


IC.+  I       lo 
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=  H 


lc+| 
DELTA  -      2.0 

TOKt   -  20.0. 

BKOlffi  -      (  £>-  l) 

ETOPL     £  ?    "  J°    1 

Input  Data.  The  data  that  is  required  is  listed  below: 
Position  Description 

0086      6  ,  total  fraction  of  delayed  neutrons. 

0106      b,  constant  of  the  reactivity  term. 

0091      (  p     -   1) 

OOlfij-      Absolute  reactivity.  Positive  value  for 
negative  insertion  of  reactivity,  f\   . 

0136      a 

0601       The  value  of  time  at  which  the  flux  is 
to       to  be  calculated,   lc. 

0621 
Output  Data.  Two  output  cards  are  punched  for  each  time  point. 
The  data  will  appear  in  the  following  form. 

Print  position.    1977   1978   1979     1980     1981   1982 
First  line.       T(C)   SUMT   SUMA-1    suma-2 
Second  line.  T(C)   NOPT 


EXPLAHAT.TON  OP  PLATE  XLV 

Plow  diagram  of  the  computer  program  for 
the  solution  of  the  reactor  kinetic  equations 
using  one  group  of  delayed  neutrons  and  when 
reactivity  is  a  function  of  time. 
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Tabulated  Data,  Analytical 

Table  5. 

,  Flux  rati 

.0  as  a  function  of 

the  value 

of  the 

decay 

constant. 

. 

Value  Used 
i'or  7.      M 
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0.6200 

0. 1)94-3 

c 
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0.10714.7 

0.01 

1.0152 

1.0272 

1.072 

1.0572 

3 

0.01 

1.0152 

1.0271 

1.0573 

1.0886 

4 

0.2900 

0.01 

1.0151 

1.0270 

1.0574 

1.0887 

5 

1.0002 

0.01 

1.0151 

1.0272 

1.0576 

1.0890 

» 

0.01 

1.0152 

1.0273 

1.0577 

1.0891 

1 

0.01216 

0.10 

Lira 

1.1736 

1.34-07 

1.9008 

2.7335 
2.7074 
2.7169 

2 

0.02958 
0 ,10747 

0.10 

1.34-01 

1.8921). 

3 

0.10 

1.1731)- 

1.3391 

1.8929 

k 

0.2900 

0.10 

1.1725 

1.3374 
I.3398 
1.3410 

I.8926 

2.7174 

5 

1.0002 

0.10 

1.1725 

I.8987 

2.7305 

« 

0.10 

1.1737 

1.901I4. 

2.7360 

Keepin's 

1  delayed  neutron  parameters. 

9   seconds 

Prompt  neutron  life 

time,  £ 

:8.0x  10 

. 

*  Here, 

all  of  the 

accepted  values  of 

the  decay  constants  were 

used. 

-:;-*  Here, 

the  value 

for  the  particular  decayed  constant 

was  the 

accepted  value 

for  the 

i-th  group 

plus  the 

experimental 

uncertainty. 
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Table  6.  Flux  ration  as  a  function  of  the  total  fraction  of 
delayed  neutrons. 


p 

:  Reac- 
:  fcivity 

1 

Plux  Ratios  at  Times  (Seo, 

.) 

;  « 

:   10 

1     l*o 

:  100   s 

1   200 

0.006 

0.01 

1.0160 

1.0292 

1.0618 

1.1191*. 

1.2112 

0.007 

0.01 

1.0131* 

1.0338 

1.0717 

1.062k 
2.021*6 

1.1391 

1.2571 

0.008 

0.01 

1.0161 

1.0295 

1.1201* 

1.2210 

0.006 

0.10 

1.1863 

1.18$ 
1.1865 

1.3736 
1.3814 

1*.0070 

12.306 

0.007 

0.10 

2.01*78 
2.0260 

i*.1021 

12.852 

0.008 

0.10 

1.3780 

1J..0128 

12.338 

0.006 

0.30 

1.81*93 

3.5377 

25.615 

1194.5 
1256.8 

6.873  x 

i°:l 
10  \ 

0.007 

0.30 

1.8558 

um3 

26 .181 

7.577  x 
7.0l6  X 

0.008 

0.30 

1.6511). 

3.51*1*2 

25.738 

1207.3 
O.3969 
O.IJ.029 

10-5 

0.006 

-0.10 

0.8622 

0.7701 

0.5977 

O.2096 

0.007 

-0.10 

O.86I4.2 

0.7722 

0.6022 

0.211*7 

0.008 

-0.10 

0.8623 

0.7703 

o.59?c 
0.1670 

C.39£l 
0.0672 

C.2C99 

0.006 

-0.50 

0.5513 

0.3795 

0.0155 

0.007 

-0.50 

0.5526 

0.3808 

0.1880 

O.0676 

0.0156 

0.008 

-0.50 

0.5517 

0.3799 

0.1873 

0.0673 

0.0155 

Hughes'  delayed  neutron  parameters.    - 

Prompt  neutron  lifetime,  _g  a  8.0  x  10--3  seconds. 
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Part  -  G 

Tabulated  Data,  Experimental 

Table  7. 

Experimental 

rod  drop  No.  A-2 

Time  after 

|   Flux 
;   Ratio 

j  Reactivity  worth 

: 

Reactivity  worth 

rod  drop 

:  of  rod  drop  (&)» 

• 

of 

rod  drop  ($), 

(See.) 

•       Hughea 

: 

Keepin 

2 

0.560 

0.545 

0.481 

j 

0.528 
0.5-79 

0.435 
0.445 

0.465 
0.471 

8 

0A55 

0.420 

0.i|.56 

10 
15 

0.3o4 

0.435 
0.425 

0.472 
O.563 

20 

0.319 

0.435 

0.468 

i 

0.259 
0.214 

0.435 
0.515 

0.463 
0.56k 

50 

0.161 

0.430 

0.466 

75 

0.121 

0.530 

0.461 

100 

0.083 

0.425 

0.458 

Reactor: 

Are 

onaut 

Critical  rod  poa 

itions: 

Shim  -  90 

Rod  dropped: 

Shim 

Table  8. 

Experimental 

rod  drop  No.  4-3. 

Time  after  : 

:   Reactivity  worth 

Reactivity  worth 

rod  drop 

Flux 
Ratio 

:   of  rod  drop  (v), 

of 

rod  drop  (£>), 

(Sec.) 

:       Hughes 

Keepin 

2 

0.521 
0.5-58 

0.560 

0.603 

i 

0.575 

0.589 

0.l|.2k 

0.550 

0.585 

8 

0.395 

0.535 

0.574 

10 

0.363 

0.535 

0.580 

15 

0.315 

0.525 

o.56k 
0.568 

£0 

0.275 

0.530 

£ 

0.219 

0.525 

0.563 

0.179 

0.525 

0.560 

50 

0.159 

0.520 

0.562 

75 

O.099 
O.069 

0.515 
0.495 

0.554 

100 

0.527 

Reactor: 

Arg 

onaut 

Critical  rod  pos 

iti 

ons: 

Shim  -  29.5 

Rod  droppi 

3d: 

Shim  and 

Fine 

Fine  -  100 
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Table  9.  Experimental  rod  drop  No.  A-I4., 


Time  after 

:                    1 

Reactivity  worth 

: 

Reactivity  worth 

rod  drop 

Flux        , 

*            D„4-»                 ' 

of  rod  drop   ((>), 

1 

of 

rod  drop   ($), 

(Sec.) 

.      Ratio      . 

Hughes 

: 

Keep in 

2 

0.956 

0.029 

0.037 

[ 

0.939 

0.033 

0.035 

0.927 

0.033 

0.03b 

8 

0.921 

0.033 

0.038 

10    ■ 

O.909 
0.838 

0.034. 

0.038 

15 

0.035 

0.038 

20 

0.872 

0.035 

0.038 

30 

O.84I 

0.03o 

0.039 

u 

0,808 

0.033 

0.039 

50 

0.796 

0.035 

O.038 

75 

0.739 
O.691 

O.O30 

0.039 
0.03B 

100 

0.035 

Reactor:     Argonaut 

Critical  rod  pos 

itions : 

Fine  -  20 

Rod  dropped: 

Fine 

Table  10. 

Experimental 

rod  drop  Ho.  A 

-5. 

Tina  after 

'      Flux 

:      Reactivity  1 

zorth     : 

Ros 

ctivity  worth 

rod  drop 

*       Ratio 

:      of  rod  droo 

tl),      : 

of 

rod  drop  ($), 

(See.) 

:                Hughsa 

Keepln 

5 

0.594 

0.541 

0.420 
0.410 

°M 

6 

o.5oo 
0.470 

0.410 

0.432 

8 

0.400 

0.433 

10 

0438 

O.40O 

0.432 

15 

0.486 

0.390 

0.1-.27 

20 

0.341 

0.282 

0.395 

0.428 

30 

0.395 

0.426 

40 

0.234 
O.I98 

0.390 

O.jj-24 

50 

0.390 

0.424 

75 

0.134 

0.395 

0.422 

100 

O.093 

0.390 

0.422 

Reactor:     j 

Irgonaut 

Critical  rod  poaiti 

ons: 

Shim  -  80 

Rod  dropped:     Shim 
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Table  11.     Experimental  rod  drop  No.  A-6. 

Time  after 

:       Reactivity  worth 

:     Reactivity  worth 

rod  drop 

.      Flint 

:       of  rod  drop   ($), 

:     of 

rod  drop   ($), 

(See.) 

.      Ratio 

[                Hughes 

• 

Keep in 

2 

0.628 

0.370 

0.393 
0.376 
0.36k 

5 

0.583 

o.3k5 

6 

0.5L7 
0.516 
O.fi-77 

0.3k0 

8 
10 

0.335 
0.3k5 

0.36k 
0.373 
0.376 

15 

0.k21 

0.3k0 

20 

0.373 

0.3k5 

0.378 

30 

0.311 

0.3k5 

0.378 

Q 

0.262 

0.3k5 

0.376 

50 

0.223 

0.3k5 

0.376 

75 

0.155 

0.3k5 

0.372 

100 

0.109 

0.3k0 

0.372 

Reactor:     Arg 

onaut 

Critical  rod  pos 

Ltions: 

Shim  -  70 

Rod  dropped: 

Shim 

Table   12.     Experimental 

rod  drop  Ho.  A-7. 

Time  after      ; 

Flux        ' 

Reactivity  worth 

:     Reactivity  worth 

rod  drop        j 

Ratio 

of  rod  drop    ($), 

:     of 

rod  drop   ($), 

(Sec.)           i 

Hughes 

• 

Keep in 

2 

0.683 

0.285 

0.306 

k 

0.618 

0.299 

0.325 

6 

0.583 
0.556 

0.2Q2 

0.318 

8 

O.283 

0.312 

10 

0.529 
0A7I 

0.281 

0.319 

15 

0.283 

0.312 

20 

O.lt-22 

0.290 

0.312 

30 

0.35k 
0.306 

0.29k 

0.316 

4-0 

0.2Q2 

0.313 

50 

0.267 

O.287 

0.307 

75 

0.188 

0.288 

0.313 
0.30k 

100 

O.ll+O 

0.282 

Reactor:     Arg 

inaut 

Critical  rod  positions: 

Shim  -  60 

Rod  dropped: 

Shim 
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Table  13.  Experimental  rod  drop  No.  A -8. 


Time  after   : 

rod  drop    • 

(Sec.)    . 


Flux 
Ratio 


0.738 
O.696 
O.661 
O.6U.3 
0.608 
0.578 


0.263 

0.201 


Reactor:     Argonaut 
Rod  dropped:     Shim 


Reactivity  worth 

of  rod  drop  (sj>), 

Hughes 


Reactivity  worth 

of  rod  drop  ($), 

Keep in 


0.221 
0.212 
0.210 
O.I99 
O.206 
0.195 
0.217 
0.209 
0.211 
0.207 
0.208 
0.217 


0.233 
0.235 
0.230 
0.219 
0.228 
0.218 
0.226 
0.226 
0.228 


0.221). 


Critical  rod  positions:  Shim  -  50 


Table  lk.     Experimental  rod  drop  No.  A-9. 


Time  after   :   p 
rod  drop    .   J*"J 
(Sec.)     I   Rati0 


0.823 
0.790 
O.768 
0.737 

o!o66 
0.632 
0.562 
O.flk 
0469 
O.382 
O.306 


Reactor:     Argonaut 
Rod  dropped:     Shim 


Reactivity  worth 

of  rod  drop    ($), 

Hughes 


0.133 
0.130 
0.125 
0.131 
0.131 
0.132 
0.131 
0.135 
0.133 
0.135 
0.135 
0.138 


Reactivity  worth 

of  rod  drop  ($), 

Keepin 


O.lkl 

o.ila 
0.136 
0.141 
o.ik5 

0.145 

0.1LL 
o.l^f 
0.114.6 
o.ilj.6 
o.lkS 


Critical  rod  positions:  Shim  -  lj.0 
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Table  IS.  Experimental  rod  drop  Ho.  A-10. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:  Reactivity  worth 
j  of  rod  drop  (£;)» 
I       Hughes 


O.892 
0.872 
0.853 
O.842 
0.822 
0.790 
0.759 


Reactor:  Argonaut 
Rod  dropped:  Shim 


0.075 
0.073 
0.071 
0.068 
0.072 
0.072 
0.07lj- 
0.074 
0.075 
0.074 
O.O76 
O.O76 


Reactivity  worth 

of  rod  drop  ($), 

Keepin 


0.081 
0.080 
0.080 
0.077 
0.080 

0.080 
0,081 
0.081 

0«08l 

0.081 
0.083 
0.083 


Critical  rod  positions:  Shim  -  30 


Table  16. 

Experimental 

rod  drop  Ho.  A- 

•11. 

Time  after 

!   Flux 
J   Ratio 

Reactivity  worth  : 

Reactivity  worth 

rod  drop 

:   of  rod  drop 

($).  : 

of 

rod  drop  ($) 

(Sec.) 

:       Hughes 

Keepin 

2 

O.968 

0.021 

C.022 

k 

0.952 
O.9E6 

0.024 

C.027 

6 

0.023 

0.027 
0.026 

8 

0.9li0 

0.023 

10 

0.931 

0.02k 

0.028 

15 

O.Qlit. 

0.026 

0.029 

20 

O.899 

0.026 

0.029 

30 

0.871 

0.028 

0.031 

40 

0.81J.9 

0.029 

0.031 

50 

0.830 

0.028 

0.031 

75 

0.778 

0.030 

0.032 

100 

0.737 

0.029 

0.032 

Reactor:  Argonaut 

Critical  rod  positi 

ons : 

Shim  -  20 

Rod  dropped 

:  Shim 
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Table  17.     Experimental 

rod  drop  No.  A- 

L2. 

Time  after 

Flux       i 
Ratio      j 

Reactivity  1 

rorth 

• 

Reactivity  worth 

rod  drop 

of  rod  drop 

(*■), 

: 

of 

rod  drop   (0), 

(See*)          i 

Hughes 

: 

Keepin 

8 
10 
15 
20 

0.988 

0.001+ 

O.OOlf 

0.983 

0.006 

0.005 

S 

§ 

100 

0.971+ 

P.001+ 

0,001+ 

0.9i4.8 

0.005 

0.005 

Reactor:     Argonaut 

Critical  rod  posi 

M 

one: 

Shim  -  10 

Rod  dropped: 

Shim 

Table   18.     Experimental 

rod  drop  Ho,  A-13. 

Time  after 

:       Reactivity  worth 

, 

Reactivity  worth 

rod  drop 
(Sec.)            j 

Flux 
Ratio 

•       of  rod  drop 
I                 Hughes 

($>. 

t 

Of    3 

rod  drop   (£), 
Keepin 

2 

0.625 
0.56L 

0.376 

taMfe 

£ 

0.375 

0405 

6 

0.516 
0.L83 

0.385 

o.l+io 

8 

0.375 

o.l+il+ 

10 
15 

0.%0 
0.386 

0.390 

0.380 

0 .400 

0.1+2$. 

0.1x26 

20 

0.31+0 

0.1+30 

30 

0.277 

o.l+oo 

0J+3p 

4.0 

0.228 

0.1+05 

0.1+31+ 

50 

0.195 

SJlQO 

0.+32 

75 

0.130 

0.505 

oJ+33 

100 

0.090 

0.1+00 

0.1+33 

Reactor:     Arg 

anaut 

Critical  roc 

positi 

nsi 

Fine  -  100 

Rod  dropped: 

Fine 
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Table  19.  Experimental  rod  drop  No.  A-l4. 


Time  after 

Reactivity 

worth 

s  Reactivity  worth 

rod  drop 

Ratio   " 

or  rod  drop 

(#), 

:  of 

rod  drop  ($) 

(See.) 

Hughes 

[ 

Keepin 

8 

0.604 

0.400 

0.435 

1 

0.553 

0.390 

0.422 

0.515 

0.475 

0.385 

0.4l0 

3 

O.39O 
0.400 

O.426 

10 

0.43? 
0.384 

0.430 
0.428 

15 

0.395 

20 
?0 

0.341 
0.275 

0.395 
0.4o5 

0.429 
0.434 

kO 

0.229 

0.4o5 

0.432 

50 

0.194 

0.400 

0.433 

75 

0.131 

0.405 

0.431 
0.426 

100 

0.091 

0.405 

Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  90 


Table  20.  Experimental  rod  drop  No.  A-15. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


Reactivity  worth 
of  rod  drop  (&), 

Hughes 


Reactivity  worth 
of  rod  drop   ($) 

Keepin 


0.683 

0.598 
0.536 
0.482 
0.438 


0.150 

0.150 

0.149 
0.150 
0.150 


0.166 

0.164 
O.163 
0.165 
O.163 


Reactor:     Argonaut 
Rod  dropped:     Fine 


Critical  rod  positions:      32.1 
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Table  21.  Experimental  rod  drop  No.  A-l6. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:   Reactivity  worth  :  Reactivity  worth 
:   of  rod  drop  ($),   :  of  rod  drop  {$) 

:       Hughes       :      Keepin 


3 

10 

0.593 

15 

20 

o.!^ 

30 

0425 

4-0 

0.375- 

50 

0.332 

75 

100 

0.218 

0.221 
0.221). 
0.22^ 
0.221 


0.2ij.3 

0.2k3 
0.2l).2 

o.aa 

0,239 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  l\.0 


Table  22.  Experimental  rod  drop  No.  A-17. 


:  txvx        s   Reactivity  worth  :  Reactivity  worth 


Time  after 

rod  drop 

(Sec.) 


Ratio 


of  rod  drop  (|i), 
Hughes 


:  of  rod  drop  ($) 
:       Keepin 


0.500 

0.39,8 
0.331). 
0.280 
0.2l<i). 


0.315 

0.316 
0.316 
0.319 
0.315 


0.314-1 

0.3W). 
0.3k0 
o.3k6 
0.3I1.I 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  50 


Table  23.  Experimental  rod  drop  No.  A-18. 
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Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


•  Keactivity  worth 

•  of  rod  drop    {£), 
:  Hughes 


Reactivity  worth 

oT  rod  drop   ($) 

Keepin 


0^36 

0.337 
0.274 
0.229 
0.1$ 


0.1+15 

0  .1+05 
0.1+05 

oJ+oo 
0 .1+00 


0^31+ 

OJ+33 
0^+3b 
0J+33 
0J+33 


Reactor:     Argonaut 
Rod  dropped :     ?ine 


Critical  rod  positions:     Fine  -  60 


Table   2i+.     Experimental  rod  drop  No.  A-19. 


Time  after 

red  drop 

(Sec.) 


Flux 
Ratio 


Reactivity  worth 

of  rod  drop   (§), 

Hughes 


Reactivity  worth 
of  rod  drop   (#) 

Keapln 


8 
10 

* 

?o 

f.0 

50 

75 
100 


0.387 
0.297 
0.237 
0.196 
0.165 


0.I+90 
0.1+80 

0.1+50 

0.L80 
0.1+75 


0.525 
0.513 
C.513 
0.510 
0.508 


Reac  tor :     Argonaut 
Rod  dropped:     Fine 


Critical  rod  positions:     iine  -  70 


Table  25.  Experimental  rod  drop  No.  A-20. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 
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Reactivity  worth  : 

of  rod  drop  (&),  : 

Hughes      j 


Reactivity  worth 

of  rod  drop  ($) 

Keepin 


0.372 


0.280 
0.225 
0.18, 
0.15 


0.520 

0.515 
0.510 
0.505 
0.500 


0.560 

0.553 
0.551 
o.&o 
53 


0.5I 


Reactor:  Argonaut 
Rod  dropped:  Pine 


Critical  rod  positions:  Fine  -  80 


Table  26.  Experimental  rod  drop  No.  A-21. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:  Reactivity  worth 
:  of  rod  drop  ($), 
:       Hughes 


Reactivity  worth 

of  rod  drop  (C>) 

Keepin 


0.365 

0.271 
0.217 
0.178 
0.lii8 


0.5I(.0 

0.535 
0.535 
0.530 
0.530 


0.575 

0.577 
0.570 
0.563 
0.562 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  90 


Table  27,  Experimental  rod  drop  No,  A-22, 


193 


Time  after 

rod  drop 

(Sec.) 


| 

8 

10 

15 


5o 
75 

100 


Flux 
Ratio 


0.362 

O.267 
0.213 


Reactivity  worth 

of  rod  drop  (;';), 

Hughe 3 


0.540 

0.550 
0.540 

0.540 

0.540 


Reactivity  worth 

of  rod  drop  (?) 

Keepin 


0.583 

0.585 
0.583 
0.584 
0.580 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  100 


Table  28.  Experimental  rod  drop  No.  A-23. 


Time  after   : 

rod  drop   j 

(Sec.)     : 


Flux 
Ratio 


Reactivity  worth 

of  rod  drop  (f?), 

Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keepin 


8 
10 

15 
20 

3° 

40 

50 

75 

100 


0.980 

0.972 
0.9&5 
0.959 
0.945 


0.007 

0.007 
0.007 
0.007 
0.007 


0.007 

0.008 
0.008 
0.008 
0.008 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  100 
Final  rod  position:   -91,1 


1* 


Table  29.  Experimental  rod  drop  No.  A-2I4.. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


Reactivity  worth  t 

of  rod  drop  ($)»   : 

Hughes       : 


Reactivity  worth 

of  rod  drop  (§) 

Keep in 


0.985 

0.938 
0.925 
0.915 
0.903 


O.OII4. 

0.016 
0.016 
0.015 
0.015 


0.016 

0.018 
0.017 
0.017 
0.017 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine 
Final  rod  position:  84.2 
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Table  30.  Experimental  rod  drop  No.  A-25. 


Time  after 

rod  drop 

(Sec.) 


Flux    s   Reactivity  worth 

Ratio   s   of  rod  dr°P  <$). 
:       Hughes 


Reactivity  worth 
of  rod  drop  ($) 

Keepin 


O.908 
0.675 
0.875 
O.8I4.3 
0.820 
0.800 


0.034 


0.033 
0.035 
0.03.' 
0.03 


0.038 

0.038 
0.039 
0.038 
0.037 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  100 
Final  rod  position:  76.5 
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Table  31.  Experimental  rod  drop  Ho.  A -26. 


Reactor:  Argonaut 
Rod  dropped!  Pine 


Time  after 

Reactivity  worth 

»  Reactivity  worth 

rod  drop 

of 

rod  drop 

(  ), 

5   Ol 

rod  drop  (t) 

(Sea.) 

Hughes 

: 

Keepln 

8 
10 

0.870 

0.050 
0.0l4.g 

0.055 
0.054 

15 

0.828 

20 

0.828 
0.792 
0.760 

0.01*9 
0.014.9 

0.0511 
0.05li- 

u 

0.050 

O.0J4.8 

0.051J. 

50 

75 
100 

0.732 

0.053 

Critical  rod  positions:  Pine  -  100 
Final  rod  positions:  73.0 


Table  32.  Experimental  rod  drop  No.  A-27. 


Time  after   :   pi     :   Reactivity  worth 
rod  drop   :   R£tl   :   of  rod  drop  ($), 
(Sec.)    :    aul°   :       Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keepln 


8 

10 

0.818 

15 

20 

0.751*. 

io 

0.702 

0.662 

50 

0.622 

75 

100 

Reactor: 

Argonaut 

Rod  dropped: 

Pine 

0.075 

0.075 
0.076 
0.076 
0.076 


0.083 


0.083 
0.085 
0.081* 
0.085 


Critical  rod  positions:  Pine 
Pinal  rod  position:  67 .1 


100 


Table  33,  Experimental  rod  drop  No,  A-28, 
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Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:  Reactivity  worth 
:  of  rod  drop  (C), 
1       Hughes 


Reactivity  worth 

of  rod  drop  (C) 

Keepin 


0.818 

0.502 
0.y6 
0.385 
0.3ijl 


0.215 

0.215 
0.215 
0.21k 
0.21k 


c;: 


0.233 
0.231 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine 
Final  rod  position:  50.0 


100 


Table  3k.  Experimental  rod  drop  No.  A-29. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:  Reactivity  worth 
:  of  rod  drop  ($), 
;       Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keepin 


6 
10 
15 
20 

& 

50 
75 

100 


0.558 

0.k56 
0.390 
0.338 
O.296 


0.252 

0.255 
0.255 
0.255 
0.25k 


0.277 
0.277 

0.277 
0.27b 
0.27k 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine 
Final  rod  position:  I+5.3 


100 
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Table  35.  Experimental  rod  drop  No.  A-30, 


Time  after 

rod  droo 

(Sec.) 


Flux 

Ratio 


Reactivity  worth 

of  rod  drop  ((; ), 

Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keep in 


O.+96 

O.4OO 
0.335 
0.286 
0.245 


0.318 

0.31)4- 
0.3l6 
0.311 
0.314 


0.314-7 

0.34? 
0.3ILO 
0.3J+0 
0.3J4.0 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine 
Final  rod  position:  39.3 


100 


Table  36.  Experimental  rod  drop  tio.  A-31. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


Reactivity  worth 
of  rod  drop  ($)# 

Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keepin 


0.^76 

0.375 
0.310 
O.263 
0.225 


0.350 

0.3U5 
0.345 
O.340 
0.345 


0.385 

0.375 
0.375 
0.37,6 
0.374 


Reactor:  Argonaut 
Rod  dropped:  Fine 


Critical  rod  positions:  Fine  -  100 
Final  rod  position:  35.2 
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Table  37.  Experimental 

rod  drop  No.  A-32. 

Time  after 

5   Flux 

!   Reactivity  worth  : 

Reactivity  worth 

rod  drop 

\      Ratio 

of  rod  drop  (£<)»   : 

of  rod  drop  (♦) 

(See,) 

Hughes       : 

Keep in 

8 
10 

O^i.0 

0.L00 

O.4.29 

15 

0.34-0 

.  0.4.00 

20 

0.34-0 

0.4.00 

0.430 

30 

0.278 

0.400 

O.430 

5-0 

0.234 
O.196 

0.399 
0.398 

0.4-25 

5o 
75 

100 

0.429 

Ileaotor:  Argonaut 

Critical  rod  poaiti 

ons:  Pine  -  100 

Rod  dropped: 

Pine 

Pinal  rod  position: 

29.2 

Table  38.  Experimental 

rod  drop  No,  A -33, 

Tine  after   : 

Flux    ' 

Reactivity  worth  : 

Reactivity  worth 

rod  drop   : 

Ratio   ! 

of  rod  drop  (§),   : 

of  rod  drop  (§) 

(Sec.)    : 

Hughes       : 

Keep in 

I 

8 

10 
15 
20 

0.370 

0.525 

O.564 

0.274- 

0.530 

0.567 

i 

0.219 

0.530 

O.563 

0.179 

0.520 

0.558 

100 

0.151 

0.520 

0.556 

Reactor:  Arg 

onaut 

Critical  rod  positi 

sns:  Pine  -  100 

Rod  dropped: 

Pine 

Pinal  rod  position: 

8.C, 

Table  39.  Experimental  rod  drop  Ho.  A-3I).. 
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Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


:  Reactivity  worth 
j  or  rod  drop  ($), 
:       Hughes 


Reactivity  worth 

of  rod  drop  ($) 

Keepin 


0.372 


0.279 
0.221). 
0.18 

0.151 


0.520 

0.515 

o.5i5 
0.505 
0.505 


0.560 

0.5A 

0.550 
0.55.0 
0.5I1.3 


Reactors  Argonaut 
nod  dropped:  Pine 


Critical  rod  positions:  Fine 
Final  rod  position:  10.3 


100 


Table  l\.0.     Experimental  rod  drop  Ho.  A-35. 


Time  after 

ioc  drop 

(Sec.) 


Reactivity  worth 
of  rod  drop  ($), 

implies 


It  activity  worth 

of  rod  drop   (#) 

Keepin 


0.383 

0.289 
0.232 
O.191 
O.161 


O.lj.95 

0495 
0.500 
0490 

04W 


0.537 

0.532 
0.530 

0.521). 

0.521 


Reactor:  Argonaut 
Kod  dropped:  Fine 


Critical  rod  positions:  Fine  -  100 
Final  rod  position:  17.0 


200 


Table  kl«  Experimental  rod  drop  No.  7. 


Time  after 

: 

PlUX     ] 

Ratio   ; 

Reactivity  worth  : 

Reactivity  worth 

rod  drop 

• 

of  rod  drop  (0),   ; 

of  rod  drop 

(  ) 

(Sec.) 

: 

Hughes       : 

Keepin 

| 

0.178 

1.75 

I.89 

e 

0.154 

1.78 

1.92 

10 

0.135 

X.79 

i.95 

15 

0.106 

1.82 

1.97 

20 

0.086 

1.83 

1.99 

Jo 

0.062 

1.85 

2.00 

75 

100 

Reactor:  TRIO* 

Critical  rod  nositi 

ons:  J  him 

-  up 

Rod  dropped* 

s 

aim 

Heg. 
Safety 

-  299 

-  up 

Table  1).2.  Experimental  rod  drop  No.  8. 


Time  after 

;  piux  | 

Ratio 

Reactivity  worth 

!  Reactivity  worth 

rod  drop 

of 

rod  drop  ($;, 

1  of  rod  drop  (♦) 

(Sec.) 

Hughes 

t      Keepin 

2 

0$30 

0.5l6 

0.560 

k 

0.515 

0.562 

5 

0.k28 

0.532 

0.578 

3 

O.396 
O.368 
O.LI4.2 

0.526 

0.574 

10 

0.52k 

0.572 

15 

0.350 

0.387 
0.386 

20 

0.390 

o*3kk 

30 

0.311 

0.375 

ko 

0.25k 

0.355 

0.387 

75 

0.212 
0.127 

O.368 

o.kio 

s$7 

100 

0.100 

o.3k5 

O.kO 

heactor:  TRICA 

hoa  dropped l  Regulating 


Critical  rod  Dosit_onsi 


-  up 

Reg.        -  308 
Safety  -  up 


Table  iv3 *     Experimental  rod  drop  No.  9, 
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Time  after 

rod  drop 

(Seo.) 

Flux 
Ratio 

f 

8 
10 

15 
20 

s 

0.258 
0.195 

0.107 
0.08Q 

0.066 

Jo 

75 

100 

0.037 
0.022 
O.Olli. 

s  Reactivity  worth 
:  of  rod  drop  (*'), 
:       Hughes 


Reactivity  worth 

of  rod  drop  (S) 

Keep in 


1.89 


1.80 
1.79 
1.77 

1.87 

1.82 
1.76 


2.02 
2.01 


I.96 
1.95 
1.92 

1.99 

1.95 
1.39 


Reactor:  TRIGA 
Rod  dropped:  Shim 


Critical  rod  positions: 


Shim 
Reg. 
Safety 


up 

299 
up 


Table  kk«  Experimental  rod  drop  Uo.  11. 


Time  after 
rod  drop 

Flux        ' 
Ratio 

Reactivity  1 
of  rod  drop 

irorth 

(    ). 

t     Res 

:     of 

ctivity  worth 
rod  droD   (*J) 

(Seo.) 

Hughes 

• 

Keep in 

2 

0.281 

1.56 

1.65 

k 

0.208 

1.75 

1.86 

0 

O.169 

1.8k 

1.99 

fl 

0.lk5 

1.85 

2.00 

10 

0.135 

1.80 

1.95 

15 

0.108 

1.78 

1.9k 

20 

0.090 
0.062 

1.76 

1.91 

30 

1.83 

1.99 

ko 

0.0l|7 

1.83 

1.99 

5o 

O.OkO 

1.77 

1.90 

75 

0.023 

1.77 

1.90 

100 

Reactor:  TRIGA 
Rod  dropped:  Shim 


Critical  rod  positions: 


Shim   -  up 
Reg.   -  311 
Safety  -  up 
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Table  4.5.  Sxperinontal  rod  drop  Ho,  13 . 


Time  after   : 

rod  drop   : 

(Sec.)    s 


Flux 


Reactivity  worth 

of  rod  drop  It), 

Hughes 


0.62$ 

0.552 
0.511 
0.4.63 
0.1+37 
0.331 
0.339 
0.27k 
0.22k 

0.1* 
0.125 

0.087 


0.365 
0.387 
0.38k 
O.4.06 
O.k.04. 

0.4.00 
0.399 
0.4-00 
0.1+13 
0.k03 

0.416 

0.1+13 


Reactivity  worth 

of  rod  drop    (( ), 

Keepln 


0.4.22 
0.423 
0.kk3 
O.W+3 
O.U33 
0.432 
O.k.27 
OA36 
0.U39 
0A3l| 
0.W 
0.4^-0 


Reactor:     TRIG-A 

Rod  dropped:     Regulating 


Critical  rod  positions: 


Shim       -  up 
Keg.       -  310 
Safety  -  up 


Table  4.6.     Experimental  rod  drop  No.   17. 


Time  after 

: 

Flux         • 
Ratio 

Reactivity  worth 

: 

Reac 

tivity  worth 

rod  drop 

: 

of 

rod  drop   (. 

). 

- 

of  r 

od  drop   (    ), 

(See.) 

: 

Hughes 

: 

Keepin 

,2 

O.4.I8 

0.85 

0.91 

i 

0.338 

0.91 

1.00 

0.300 

0.90 

O.98 

3 

0.258 

6.9k 

1.03 

10 

0.236 

O.94. 

1.03 

15 

0.196 

0.164. 

0.93 
0.9k 

1.01 

20 

1.01 

30 

0.127 

0.9k 

1.01 

4.0 

0.101 

0.93 

1.00 

50 

0.082 

0.92 

0.98 

75 

0.050 

O.QO 

0.96 

100 

0.033 

O.89 

0.92 

Reactor:     TRIGA 

C; 

■itical  rod 

pos 

iti 

ons: 

Shim       -  380 

Rod  dropped! 

Shim 

Reg.       -  577 

Safety  -  up 

ao3 


Table  1).7.     Experimental  rod  drop  So.  20. 


Time  after 

1      Tluy         s 

Reactivity  worth 

j     Reactivity  worth 

rod  droo 

5      Ratio      ! 

of 

rod  drop 

(*). 

•     of 

rod  drop   (9), 

(Sec.) 

Buphes 

: 

Kaeplrs 

2 

0.289 

1.50 

1.6l 

1 

8 
10 

0.218 

1.65 

1.78 

0.187 

1.65 

1.80 

15 

0.113 

1.70 
l.§7 

1.80 

20 

0.0Q5 
O.069 

1.80 

30 

1.66 

1.83 

ko 

CO  53 
o.oJ+i 

1.70 

1.82 

50 

1.73 

1.7k 

1.81 

75 

0.023 

1.83 

Reactor:     TRIOA 

Rod  dropped:     Regulating 


Critical  rod  positions: 


3  aim 
i'^eg. 
Safety 


up 


Table  1).8.  Experimental  rod  drop  Ko.  21. 


Time  after 

, 

Reactivity  worth      : 

heac 

tivity  worth 

rod  drop 

Flux 

of 

rod  drop   [\ 

).     1 

of  1 

od  drop 

(♦). 

(See.) 

.       Ratio 

Hughe  s 

Keepin 

2 

0438 

0.79 

0.81). 

i 

O.36I 
0.315 

0.83 
0.84 

0.90 
0.92 

8 

0.285 

0.63 

0.91 

10 

15 

0.229 

0.80 

0.85 

20 

0.187 

0.82 

0.89 

3° 

0.114.2 

O.83 
0.6k 

0.90 

to 

0.112 

0.91 

50 

0.088 

0.85 

0.91 

75 

0.055 
0.034 

0.83 
0.64 

0.Q0 
O.89 

ICO 

Reactor:     TRIGA 

Critical  rod 

positi 

oris : 

Shim 

-  up 

Rod  cropped: 

Regulating 

Seg. 
Safety 

-438 
-  up 

20lf 

Table  1).9.  Experimental 

rod  drop  No.  16. 

Time  after 

TlMK 

!   Reactivity  worth 

1 

Reactivity  worth 

rod  drop 

Ratio 

:   of  rod  drop  ( & ) , 

• 

of 

rod  drop 

(  ). 

(See.) 

Hughes 

• 

Keep in 

2 

0.356 

1.10 

1.18 

J 

0.283 

1.17 

1.27 

0.2l|2 

1.20 

1.31 
I.36 

8 

0.207 

1.25 

10 

15 

0.150 

1.25 

1.35 
1.3& 

20 

0.126 

1.25 

30 

0.097 
0.074. 

1.21 

I.24 

1.36 

i3fc 

75 

O.OltO 

1.2& 

1.36 

100 

Reactor- 1  T2IGA 

Critical  rod  poa 

iti 

ona 

3'alia 

-  466 

Rod  dropped: 

Shim 

Reg. 
Safety 

-klh 
-  up 

Table  $0.     Experimental 

rod  drop  No.  19. 

Tine  after 

j  'lux 

Reactivity  worth 

1 

Reactivity  worth 

rod  drop   s 

Ratio 

of  rod  drop  ($), 

• 

of 

rod  drop 

($), 

(Sec.)'    i 

Hughes 

: 

Keep in 

8 
10 

15 

0.16% 

1.06 

I.09 
1.06 

20 

0.157 

0.99 
O.96 

30 

0.123 

1.03 
l.Olj. 

M 

O.O55 
0.088 

O.96 

5o 

0.92 

1.00 

75 

C.051 

0.94 
O.96 

1.02 

10C 

O.036 

O.98 

Reactor:  TRI 

GA 

Critical  rod  pes 

ltions: 

Shim 

-  up 

Rod  dropped: 

Regulating 

Reg. 

-  519 

Safety 

-  up 

205 


Table  51, 

Experimental  rod  drop  No,  22. 

Time  after 

!    rtm     s 

Reactivity  worth 

« 

Reactivity  worth 

rod  drop 

•                 iLL*-               • 

of  rod  drop    (£ ), 

1 

of 

rod  droo    (?), 

(See.) 

Hughes 

: 

Keepin' 

2 

oJOi-3 

0.367 

0.78 

0.82 

£ 

0.81 

0.32 

6 

0.318 

O.83 

O.91 

8 

O.290 
0.261 

0.82 

O.89 

10 

O.83 

0.8k 

O.91 

15 

0.213 

O.92 

20 

0.183 

0.8k 

0.91 

30 

O.lkO 

0.84 

0.91 

5o 

0.109 
0.086 

0.65 
0.82 

0.92 
O.93 

75 

0.053 

0.86 

0.91 

100 

0.034 

0.8k 

0.90 

Reactor j     TRIGA 

Critical  rod  positi 

ons: 

Shim       -  up 

Rod  dropped 

:     Regulating 

Reg.       -  k38 
Safety  -  up 

Table  52. 

Experimental 

rod  drop  No.  2k, 

Time  after 

s       Flux 

j      Reactivity  worth 

• 

Reactivity  worth 

rod  drop 

:       Ratio 

:       of  rod  drop   ($), 

5 

of 

rod  drop   (?>)» 

(Sec.) 

:                Hughes 

1 

Keepin 

1 

O.526 

0.55 
0.59 

0.60 
O.63 

O.kOO 

0.59 

0.65 

8 

0.357 

0.6l 

0.67 
0.66 

10 

0.330 

0.6l 

15 

0.276 

O.62 

0.67 

20 

0.237 

0.62 

0.68 

30 

0.179 

0.65 

0.70 

M 

0.lk2 

0.66 

0.71 

50 

0.116 

0.66 

0.71 

75 

0.071 

0.66 

0.71 

100 

Reactor i     TRIGA 

Critical  rod  positi 

ons: 

Shim       -  296 

Rod  dropped 

i     Shim 

Reg,       -  down 
Safety  -  3/k 
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Table  53. 

Experimental  rod  drop  Ho.  26. 

Time  after 

I       Flux 

1       Ratio      . 

Reactivity  worth 

t 

Reactivity  worth 

rod  drop 

of  rod  drop   ('•  ), 

of 

rod  drop   ($), 

(Seo,) 

Hughes 

: 

Ksepin 

I 

0.339 

1.23 

1.26 

6 
10 

0.163 

I.I4.O 

1.51 

15 

0.133 

I.I42 

1.514 

20 

30 

0.082 

1.1*3 

1.52 

M 

O.063 

143 

1.514 
1.54 

5o 

0.050 

IM 

75 

0.029 

1.I42 

1.53 

100 

Reactor:     ' 

rRIGA 

Critical  rod  positions: 

Shim       -  500 

Rod  dropped;     Shim 

Reg.        -  Clown 

Safety  -  l/2 

Table  5I4. 

Experimental  rod  drop  No,  29»1. 

Time  after 

1       Flux         ! 
J      Ratio      ! 

Reactivity  worth 

t 

Reactivity  worth 

rod  drop 

of  rod  drop   ($), 

t 

of 

rod  drop   (v), 

(Sec.) 

Hughes 

: 

Keepin 

{ 

0.283 

1.56 

1.65 

6 

8 

10 

15 

0.1014. 

1.85 

1.98 

20 

0.085 

1.85 

2.00 

30 

O.063 

1.85 

2.00 

U 

0.0l48 

1.87 
1.86 

2.01 

50 

O.038 

2.00 

75 

0.022 

1.65 
I.8I4 

1.99 

100 

0.013 

1.97 

Reactor:     TRIOA 

Critical  rod  positions: 

Shim       -  338 

Rod  droppei 

i:     Regulating 

Reg.        -  805 
Safety  -  down 
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Table  55. 

Experimental  rod  drop  No.   29-2. 

Time  after 

1       Flux        ! 
J       Ratio      J 

Reactivity  worth 

• 

Reactivity  worth 

rod  drop 

of  rod  drop   ($), 

J 

of 

rod  drop   ((•;), 

(See.) 

Hughes 

• 

Keep in 

2 

1 

0.283 

1.56 

1.65 

0.171 

1.8k 

2.00 

8 

0.1M 

1.85 

2.00 

10 

0.131 

1.85 

2.00 

15 

0.10k 

1.85 

1.98 

20 

0.085 

1.85 

2.00 

?° 

0.063 

1.85 

2.00 

Eg 

0.0k8 

1.87 
1.86 

2.01 

52 

0.036 

2.00 

75 

0.022 

1.85 
1.8k 

1.99 

100 

0.013 

1.97 

Reactor: 

EUflA 

Critical  rod  positions: 

Shim       -  338 

Rod  droppe 

d:     Regulating 

Reg.       -  805 
Safety  -  down 

Table  56. 

Experimental  rod  drop  Ho.  37-1. 

Time  after 
rod  drop 

!      Flux        i 
•       Ratio 

Reactivity  worth 
of  rod  drop   ($), 

i 
t 

Reactivity  worth 
of  rod  drop   (0)# 

(Sec.) 

•                        • 

Hughes 

: 

iicepin 

2 

0.326 

1.27 

1.31 

1.&5 

1 

0.257 

i.3k 

0.218 

1.37 

142 

1.50 

8 

0.192 

1.38 

10 

15 

0.130 

1.1)8 

1.60 

20 

?° 

0.076 

1.53 
1.56 
1.56 

1.66 

i 

0.058 
0.0k5 

1.67 

1.68 

75 

0.025 

1.57 

1.56 

I.65 

100 

0.016 

1.6? 

Reactor:     TRIGA 

Critical  rod  positions: 

Shim       -  k23 
Heg.        -  598 

Rod  dropped:     Regulating 

Safety  -  down 

2l'c. 


Table  57.  Experimental  rod  drop  Ho.  37-2. 


Time  after 

rod  drop 

(Sec.) 


Flux 
Ratio 


Reactivity  worth  : 

of  rod  drop  (<•),   : 

Hughes      : 


Reactivity  worth 

of  rod  drop  (*)» 

Keepin 


0.3: 

0.2* 


1»23 

1.29 


1.32 

1.40 


0.16I). 

0.10l|. 

0.077 

O.OpS 
0.04$ 

0.020 

0.012 


1.39 

1.52 
1.50 
1.56 
1.57 
1.55 
1.55 


1.52 

1.65 

iM 
1.67 
1.69 
1.66 

1.63 


Reactor:  IBIGA  Critical  rod  positions:  Shim   -  423 

Rod  dropped:  Regulating  Keg.   -  596 

Safety  -  down 


Table  58.     Experimental  rod  drop  No.  39. 


Time  aftsr 

Flux    | 
Ratio   : 

Reactivity  worth 

j  Reactivity  worth 

rod  drop 

of 

rod  drop  [v)t 

:  of 

rod  drop  ($), 

(Sea.)    1 

Hughes 

: 

Keepin 

,2 

O.360 

1.09 

1.17 

fc 

0.286 

1.15 

1.18 

6 

0.245 

1.17 

1.26 

8 

0.216 

1.19 

1.27 

10 

0.197 

I.16 

1.30 

15 

0.154 

1.20 

1.24 

20 

0.131 

1.20 

1.30 

30 

O.O98 

1.19 

1.30 

4-0 

0.078 
O.062 

1.17 

I.26 

50 

1.17 

1.25 

75 

100 

Reactor:     TRIGA  Critical  rod  positions:     Shim  -  495 

Rod  dropped:     Regulating  Reg.  -  520 

Safety  -  down 
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Table  59.  Experimental  rod  drop  Ho.  27. 


Time  after 

t 

Flux 
Ratio 

:      Reactivity  worth 

•     Reactivity  worth 

rod  drop 

t 

!       of 

rod  drop 

c:. 

of 

rod  drop   (:*'■)» 

(Sec.) 

: 

Hughes 

Xeepin 

0.282 

1.56 

1.65 

k 

0.231 

1.57 

1.66 

0 

0.195 

1.58 

1.70 

8 

0.170 

1.59 

1.71 

10 

0.114.5 

1.65 

1.80 

15 

0.107 

1.80 

!•% 

20 

O.O89 
O.066 
O.OSO 
O.OkO 

1.78 

1.79 
1.80 

1.93 
i.9k 
1.94 

50 

1.79 

1.92 

75 

0.021 

1.53 
1.8k 

1.99 

100 

0.012 

1.99 

Reactor:     IRICfA 

Critical  rod 

positi 

MM 

Shim  -  600 

Rod  dropped: 

Shim 

Peg.  -  down 

Table  60.  Experimental  rod  drop  No.  28. 


Time  after 

Flux 
Ratio 

:       Reactivity  \ 

vorth 

1 

reactivity  worth 

rod  drop 

(Sec.) 

: 

:       of  rod  drop 
1                 iiugUes 

(•   )» 

t 
t 

of 

rod  drop    (v)# 
iCeepin 

t 

0.292 

1A9 

1.57 

6 

e 

0.150 

1.31 

1.98 

x2 

0.136 

1.32 

51 

1.90 

15 

0.107 

1.S3 

20 

0.085 

1.33 

1.97 

30 

0.06k 
o.ok9 

1.86 

1.38 

1.98 
1.96 

50 
75 

0,037 
0.022 

1  ?7 

1.99 
1.96 

100 

Reactor:     TRIGA 

Critioal  rod 

DOSitJ 

ons: 

Shim       -  693 

Rod  dropped: 

c 

him 

Reg,       -  down 
Safety  -  l/k 
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Table  6l.     * 

Ixperimente.l  rod  drop  No.  36. 

Time  after 

*       Flux         " 

Reactivity  worth     t 

Reactivity  worth 

rod  droo 

Ratio 

of  rod  drop   (^),     : 

of 

rod  drop   ($), 

(See.) 

Hughes                : 

Keep in 

2 

0.523 

0J& 

0.55 

0.59 

1 

-..- 

0.62 

04(4 

0.58 

0.61+ 

e 

0.372 

0.58 

0.63 

10 

0.372 

0.59 

0.60 

0.63 

is 

0.292 

0.62 

20 

0.256 

O.67 

0.66 

P 

O.I89 
0.4.8 

o.6i 

0.66 

to 

0.63 
0.64 

0.67 

5o 

0.122 

0.67 
0.66 

75 

0.077 

0.62 

100 

0.051 

0.62 

0.61+ 

Reactor:     THIGA 

Critical  rod  positions: 

Shim       -  389 

Rod  dropped: 

Regulating 

heg.       -  389 
Safety  -  down 

Table  62.     I 

sperlmental  rod  drop  No.  A-l. 

Time  aftor 

•       Flux 

Reactivity  worth     : 

Reactivity  worth 

rod  drop 

|       Ratio       | 

of  rod  drop   (v  ),      : 

of 

rod  drop   (    ), 

(See.) 

Hughes                : 

Keep in 

\ 

0.5J+9 
0.510 
0.(1.69 

O.505 
0465 

0.514 
0.501 

O.IjOO 

6 

0465 

8 

04-37 

0455 

0490 

10 

o.lj.02 

0460 

0497 
0496 

IS 

0.3i)-7 

o46o 

20 

0.309 

0455 

0490 

30 

0.250 

0455 

O.lj.87 
0476 

ho 

0.208 

o45o 

75 
100 

o.nk 

o45o 

0484 

Reao  tor :     Argonaut 

Critical  rod  positions 

!     Shim  -  100 

Rod  dropped: 

Shim 
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An  analytical  study  was  made  of  the  reactor  kinetic  equa- 
tions for  a  bare  homogeneous  thermal  reactor  assuming  step 
changes  of  reactivity.  The  study  considered  two  sets  of  delayed 
neutron  parameters,  Hughes'  and  Keepln's,  and  the  results  of  the 
study  are  given  in  graphical  form.  Experimental  measurements  of 
the  reactivity  worth  of  the  control  rods  in  the  Torry  Pines 
TRIGA  reactor  were  made.  The  reactivity  values  obtained  with 
positive  period  measurements  were  compared  to  those  obtained  by 
means  of  rod  drop  measurements.  Agreement  of  the  two  methods 
indicated  that  the  kinetics  theory  was  valid  when  the  reactor 
was  operating  at  low  power  levels  and  reactivity  changes  were 
not  extremely  large  (  f>   <  $2.00). 

A  parameter  analysis  of  the  kinetic  equations  was  made. 
The  reactor  parameters,  prompt  neutron  lifetime,  total  fraction 
of  delayed  neutrons,  and  decay  constants  of  the  delayed  neu- 
trons, were  evaluated  over  a  relatively  wide  range  of  values. 
The  importance  of  each  parameter  in  the  kinetic  equations  is 
given. 

An  approximate  solution  to  the  kinetic  equations  was  ob- 
tained using  three  groups  of  delayed  neutrons  instead  of  the 
usual  six  groups.  The  delayed  neutron  parameters  were  empiri- 
cally chosen  to  fit  the  computed  data  obtained  with  six  delayed 
neutron  groups . 

Kinetics  theory  for  a  finite  reactivity  insertion  rate  was 
developed.  The  kinetics  model  included  the  assumption  that  the 
delayed  neutrons  could  be  approximated  by  one  group  of  delayed 


neutrons.     The  reactivity  insertion  rate  was  approximated  by  a 

-hfrn 

function,  of  the  form  A(l  -  e  "  ),  to  fit  the  reactivity  inser- 
tion rate  of  the  TRIGA  reactor.  The  finite  insertion  rate  equa- 
tions predicted  the  neutron  flux  to  a  greater  degree  of  accuracy 
for  short  times  after  the  reactivity  change,  than  the  neutron 
flux  as  predicted  by  the  usual  six  group  kinetic  equations  with 
step  input  reactivities. 

Three  IBM-650  programs  are  given,  solution  of  the  kinetic 
equations  for  six  groups  of  delayed  neutrons,  solution  of  the 
kinetic  equations  for  three  groups  of  delayed  neutrons,  and 
solution  of  the  one  delayed  neutron  group  kinetic  equations  when 
reactivity  is  a  function  of  time. 


